Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  453 1707 (N20E64 L=131)  1526    292        3 ( 429/1662- 498/1756) -41.2 E70-E59 N16-N23
      2  701 1004 (S16E18 L=177)   497    131        0 ( 639/ 969- 755/1044)  50.2 E21-E15 S19-S12
      3 1836  832 (N09W60 L=255)  1313    203        5 (1800/ 769-1865/ 896) -33.7 W55-W66 N06-N12
      4  829  827 (S22E06 L=189)   238     53        0 ( 803/ 802- 858/ 858) -14.0 E07-E03 S23-S21
      5 1427  773 (S08W30 L=225)   624    114        2 (1408/ 712-1448/ 835) -66.8 W27-W32 S11-S04
      6 1013  754 (S21W07 L=202)   351    149        0 ( 953/ 726-1093/ 795)   5.7 W02-W12 S22-S19
      7 1376  696 (S13W30 L=225)   236     71        2 (1343/ 668-1399/ 727) -14.0 W27-W31 S14-S12
      8  609  631 (S40E16 L=179)   771    205        1 ( 513/ 611- 707/ 654)  14.0 E23-E07 S43-S37
      9 1354  462 (S27W39 L=234)  1329    469        1 (1132/ 430-1523/ 514)  19.7 W23-W51 S32-S20
     10 1650  419 (S18W65 L=260)   497    160        1 (1607/ 391-1681/ 434)  12.5 W59-W68 S19-S17
     11 1226  253 (S43W47 L=242)   455     84        0 (1190/ 234-1264/ 269)   9.5 W42-W49 S44-S42
  

    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