Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  883 1756 (N58W03 L= 62)   379     54        1 ( 854/1733- 915/1773)  -9.5 W00-W06 N56-N59
      2  687 1721 (N59E21 L= 38)   392     55        1 ( 657/1701- 710/1740)  18.4 E25-E17 N56-N60
      3 1065 1439 (N31W09 L= 68)   328    123        0 (1018/1392-1091/1464) -63.4 W06-W11 N27-N33
      4  852 1371 (N30E06 L= 53)   376    127        0 ( 806/1339- 914/1414) -35.5 E09-E02 N27-N33
      5 1558 1308 (N15W40 L= 99)   259     55        1 (1538/1277-1585/1329) -90.0 W38-W41 N12-N15
      6 1567 1242 (N11W38 L= 97)   325     63        0 (1533/1217-1598/1264) -26.6 W36-W40 N08-N12
      7  414 1127 (N19E40 L= 19)   313    104        1 ( 400/1086- 432/1174) -80.5 E41-E38 N16-N22
      8 1903  994 (S11W67 L=126)   565     50        0 (1890/ 964-1918/1022)  21.8 W64-W70 S13-S09
      9 1563  960 (S06W33 L= 92)   223     72        0 (1552/ 922-1574/1001) -90.0 W32-W33 S08-S03
     10  185  945 (N09E64 L=355)   471     61        0 ( 168/ 915- 204/ 977)  53.1 E66-E61 N07-N11
     11 1229  751 (S13W09 L= 68)   345    112        0 (1189/ 727-1282/ 786)  18.4 W06-W12 S14-S11
     12 1759  616 (S33W52 L=111)  1996    272        0 (1668/ 526-1865/ 688)  14.0 W43-W68 S37-S29
     13  394  581 (S14E51 L=  8)   831    187        0 ( 348/ 504- 420/ 665)  42.5 E59-E46 S19-S08
     14  247  603 (S12E68 L=351)   612     58        0 ( 224/ 579- 274/ 619) -26.6 E71-E63 S14-S11
     15 1281  365 (S42W09 L= 67)   480     79        0 (1241/ 345-1326/ 383)   8.1 W05-W12 S43-S40
     16 1091  270 (S46E09 L= 50)  1134    132        1 (1034/ 234-1167/ 298)   4.8 E14-E02 S49-S43
     17  926  203 (S50E27 L= 31)  1349    139        1 ( 886/ 151- 967/ 258)  40.6 E36-E20 S56-S44
  

    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