Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  788 1346 (N12E21 L=103)   308     89        0 ( 760/1320- 833/1381)  -9.5 E23-E17 N10-N13
      2  679 1143 (S02E23 L=101)   473    185        0 ( 627/1060- 742/1220) -31.0 E27-E17 S05-N01
      3 1891 1105 (N22W62 L=187)  1246    218        3 (1842/1069-1919/1153)  40.6 W56-W66 N19-N25
      4  253 1105 (S12E55 L= 69)   569     79        1 ( 231/1062- 265/1151) -68.2 E57-E53 S14-S09
      5  921 1039 (S03E06 L=118)   284     97        1 ( 898/ 989- 942/1080)  90.0 E07-E05 S06-N00
      6 1853  993 (N14W60 L=184)   395     85        0 (1827/ 975-1877/1009)  18.4 W56-W62 N13-N15
      7  680  898 (S16E18 L=106)   897    285        0 ( 548/ 812- 773/1003)  51.3 E26-E14 S23-S08
      8 1392  879 (S03W25 L=149)   233     47        1 (1378/ 849-1406/ 917) -53.1 W23-W26 S04-N00
      9 1519  692 (S11W39 L=163)   276     61        0 (1482/ 673-1552/ 716)  45.0 W36-W40 S12-S08
     10 1308  581 (S22W28 L=153)   234     68        0 (1268/ 562-1343/ 608)  53.1 W26-W29 S23-S19
     11 1709  533 (S16W62 L=186)   589     78        0 (1677/ 505-1741/ 558)  38.7 W59-W64 S17-S13
     12 1216  368 (S38W32 L=156)  1098    264        0 (1085/ 347-1341/ 386)  22.4 W22-W40 S41-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