Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1246 1852 (N65W41 L=348)  1350    152        2 (1178/1829-1309/1886)   6.0 W28-W67 N62-N67
      2 1536 1574 (N37W47 L=354)   340     48        1 (1507/1558-1567/1587)   0.0 W43-W50 N36-N38
      3  926 1378 (N27E06 L=301)  1060    388        2 ( 867/1292-1031/1470) -57.5 E09-W01 N21-N33
      4 1535 1382 (N24W38 L=345)   363     69        0 (1512/1350-1555/1416) -76.0 W37-W39 N22-N26
      5 1449 1225 (N15W29 L=336)   268     75        0 (1423/1192-1479/1248)  31.0 W26-W31 N12-N16
      6  218 1101 (N09E61 L=246)   715     98        2 ( 199/1045- 237/1151) -66.8 E63-E58 N05-N12
      7 1158 1137 (N11W09 L=316)   194     28        1 (1131/1118-1185/1151) -18.4 W07-W10 N09-N11
      8 1747  856 (S10W51 L=358)   628     18        0 (1723/ 828-1772/ 883)  71.6 W48-W53 S11-S08
      9  767  814 (S08E17 L=290)   699    172        1 ( 694/ 777- 826/ 843)  -7.1 E21-E13 S10-S06
     10 1258  656 (S20W14 L=321)   241     73        2 (1232/ 621-1279/ 695)  68.2 W12-W15 S22-S17
     11  318  651 (S19E56 L=251)   573     72        1 ( 289/ 619- 355/ 676)  12.5 E59-E50 S21-S17
     12 1532  589 (S26W35 L=342)   403     82        1 (1507/ 565-1568/ 617) -31.0 W33-W39 S28-S24
     13 1166  151 (S67W16 L=323)   690     48        0 (1152/ 120-1177/ 182) -49.8 W12-W25 S74-S61
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y