Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  752 1802 (N63W09 L= 67)   359     30        1 ( 735/1780- 768/1824)  80.5 W05-W11 N60-N66
      2  656 1751 (N62E06 L= 52)   439     50        1 ( 627/1726- 684/1772) -14.0 E09-E01 N59-N64
      3  461 1667 (N62E34 L= 24)  1579    208        3 ( 392/1607- 556/1708)   0.0 E47-E19 N57-N65
      4 1531 1268 (N05W36 L= 94)   312    119        0 (1479/1247-1584/1296) -33.7 W33-W39 N02-N07
      5 1791 1235 (S04W55 L=113)  1158    152        3 (1739/1179-1846/1268)   9.5 W49-W61 S06-S02
      6 1854 1148 (S11W59 L=117)   472     90        0 (1828/1117-1891/1184) -39.8 W56-W63 S13-S08
      7  330 1111 (N25E44 L= 14)   604     75        1 ( 300/1072- 362/1159)  59.0 E46-E40 N23-N28
      8  849 1122 (N14E07 L= 50)   287     97        2 ( 792/1103- 889/1148)  -9.5 E10-E04 N12-N15
      9 1116 1119 (N08W07 L= 65)   192     57        1 (1085/1102-1146/1142) -18.4 W05-W08 N06-N08
     10 1112 1070 (N05W06 L= 64)   235     13        0 (1088/1049-1131/1091) -45.0 W04-W07 N03-N06
     11  901  821 (S04E12 L= 46)   988    299        1 ( 788/ 744-1004/ 930)  16.4 E20-E03 S05-N00
     12 1690  575 (S41W36 L= 94)   574     76        0 (1658/ 544-1725/ 616)  31.0 W32-W39 S43-S39
     13  472  516 (S11E49 L=  9)   355     68        0 ( 458/ 479- 494/ 543)  31.0 E51-E46 S13-S09
     14  555  511 (S13E43 L= 15)   255     61        0 ( 541/ 484- 574/ 542)  76.0 E44-E41 S15-S11
     15  339  442 (S14E68 L=350)   661     76        0 ( 324/ 403- 355/ 485)  45.0 E71-E63 S16-S10
     16 1409  393 (S45W08 L= 66)   301     55        1 (1387/ 374-1436/ 411)   0.0 W06-W10 S46-S43
     17 1152  298 (S44E15 L= 43)   563    124        1 (1088/ 280-1218/ 323) -15.3 E19-E08 S44-S41
     18  927  226 (S42E35 L= 23)   477     81        0 ( 890/ 201- 969/ 251)   6.3 E39-E30 S43-S40
  

    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