Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1619 1575 (N18W58 L=224)  3298    575        7 (1471/1538-1780/1611) -22.8 W46-W78 N11-N24
      2 1496 1354 (N11W37 L=204)   342    107        0 (1480/1297-1509/1404)  63.4 W35-W40 N07-N13
      3 1898 1203 (S11W67 L=234)  1747    301        7 (1870/1137-1936/1264)  60.3 W61-W75 S14-S07
      4  588 1170 (N26E23 L=144)   843    232        1 ( 495/1074- 655/1244) -60.3 E28-E19 N18-N32
      5 1833  959 (S22W51 L=218)  2995    357        2 (1778/ 794-1862/1122)  79.2 W46-W54 S31-S10
      6  137  932 (N21E67 L= 99)   703    114        0 ( 127/ 874- 147/ 990)  60.3 E69-E64 N17-N24
      7 1686  782 (S27W35 L=201)  1150    282        1 (1661/ 677-1734/ 892)  77.0 W32-W39 S35-S20
      8 1365  747 (S18W12 L=178)   248     87        0 (1322/ 712-1396/ 782) -71.6 W10-W13 S20-S14
      9 1062  630 (S17E09 L=158)   880    265        1 ( 975/ 516-1138/ 745)  33.7 E17-E01 S21-S10
     10  568  607 (S06E39 L=127)  1118    181        2 ( 517/ 531- 635/ 686) -84.8 E41-E36 S11-N00
     11 1196  451 (S32E07 L=160)   523    102        2 (1162/ 416-1254/ 488)  14.0 E09-E01 S33-S30
     12 1003  343 (S34E23 L=143)  1291    231        2 ( 928/ 265-1089/ 410)  16.4 E32-E15 S37-S29
  

    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