Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  804 1618 (N33E18 L=294)   289     39        1 ( 787/1592- 828/1644) -71.6 E19-E16 N31-N34
      2 1183 1525 (N27W09 L=320)   338     91        0 (1151/1497-1226/1562) -38.7 W06-W11 N25-N29
      3 1431 1321 (N15W25 L=336)   721    302        0 (1296/1253-1558/1361) -12.5 W15-W33 N11-N17
      4 1756 1270 (N14W50 L=  1)  1474    276        1 (1730/1158-1812/1360) -90.0 W47-W57 N07-N20
      5 1394 1107 (N02W22 L=333)   345     95        0 (1373/1067-1414/1160) -78.7 W20-W23 N00-N05
      6 1921 1027 (N01W67 L= 19)  1261    100        0 (1871/ 994-1979/1053)   6.0 W60-W79 S01-N03
      7 1817  937 (S05W55 L=  7)   429     78        0 (1785/ 897-1853/ 966)  38.7 W52-W59 S07-S03
      8 1829  800 (S13W59 L= 11)  2444    250        8 (1767/ 698-1902/ 919)  90.0 W54-W68 S19-S06
      9 1731  830 (S12W49 L=  0)   347     80        0 (1718/ 790-1747/ 865) -76.0 W47-W50 S14-S10
     10  803  674 (S26E13 L=298)  1128    341        0 ( 732/ 566- 932/ 762)  45.0 E18-E04 S32-S19
     11  396  717 (S23E43 L=268)   398     77        1 ( 372/ 681- 427/ 757) -78.7 E45-E40 S25-S20
     12 1750  685 (S21W55 L=  6)   340     35        0 (1727/ 668-1776/ 703)  -8.1 W51-W58 S21-S19
     13 1513  614 (S27W36 L=347)  1025    175        2 (1417/ 589-1599/ 646)  12.1 W28-W42 S28-S24
     14  864  551 (S34E10 L=302)   204     35        0 ( 845/ 534- 887/ 575) -45.0 E11-E08 S34-S32
     15 1748  459 (S34W70 L= 21)   767     84        0 (1711/ 438-1792/ 477)   0.0 W62-W78 S35-S32
     16 1372  391 (S43W32 L=343)  2465    445        3 (1143/ 367-1573/ 432)   7.7 W12-W50 S45-S38
     17  638  361 (S48E33 L=278)   653     82        2 ( 604/ 337- 687/ 390) -26.6 E37-E28 S49-S45
     18  503  325 (S51E53 L=259)   709     96        1 ( 453/ 308- 541/ 360)   4.1 E61-E45 S52-S48
  

    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