Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1083 1476 (N27W13 L=212)   597    167        1 (1009/1426-1154/1516) -53.1 W08-W16 N22-N30
      2  983 1420 (N25W05 L=204)   193     85        0 ( 939/1397-1021/1436) -31.0 W02-W07 N23-N26
      3 1543 1370 (N11W39 L=238)   737    144        0 (1481/1324-1611/1403) -45.0 W35-W42 N07-N14
      4  842 1215 (N16E07 L=192)   427    187        0 ( 768/1170- 919/1279)  20.0 E12-E01 N14-N18
      5  135 1133 (N23E67 L=132)   501     82        1 ( 103/1111- 166/1148) -10.3 E72-E61 N21-N24
      6 1870 1098 (S10W61 L=260)  1355    120        5 (1830/1049-1921/1148)  51.3 W55-W67 S13-S07
      7  197 1048 (N17E58 L=141)  2215    159        3 ( 146/ 979- 276/1116)  45.0 E64-E49 N13-N21
      8 1270 1021 (S03W14 L=213)   287    111        0 (1217/ 990-1309/1039) -31.0 W10-W16 S05-S01
      9  260  974 (N12E52 L=147)   459     83        2 ( 242/ 938- 280/1012)  53.1 E53-E49 N10-N14
     10 1172  899 (S08W06 L=205)   497    140        0 (1146/ 848-1215/ 968) -81.9 W04-W07 S11-S04
     11 1320  872 (S13W14 L=213)   355    161        1 (1248/ 834-1377/ 930)  18.4 W09-W18 S13-S10
     12  385  577 (S13E52 L=147)  2043    428        3 ( 330/ 444- 421/ 788)  46.5 E65-E45 S20-N00
     13 1155  756 (S16W03 L=202)   225    103        0 (1112/ 721-1187/ 782) -53.1 W00-W04 S18-S13
     14 1685  603 (S37W41 L=240)  3214    645        2 (1428/ 434-1910/ 780)  14.6 W16-W68 S44-S30
     15  606  648 (S12E33 L=166)  1291    327        3 ( 480/ 589- 718/ 745)  14.7 E42-E23 S14-S08
     16  181  644 (S06E73 L=126)  1102    111        0 ( 151/ 597- 208/ 696) -71.6 E76-E70 S08-S02
     17  595  516 (S20E38 L=161)   533    135        1 ( 566/ 476- 621/ 555)  68.2 E41-E35 S22-S17
     18  534  410 (S26E49 L=150)   976    161        1 ( 458/ 382- 625/ 451)   0.0 E56-E39 S27-S23
     19  923  364 (S37E23 L=176)   252     20        0 ( 905/ 345- 941/ 385) -90.0 E24-E21 S38-S34
  

    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