Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  837 1697 (N30E30 L= 16)  1968    445        4 ( 670/1578- 982/1800) -14.8 E49-E15 N24-N35
      2 1265 1753 (N46E02 L= 44)   579     91        1 (1227/1717-1295/1788)  69.4 E04-W00 N42-N50
      3  462 1685 (N23E59 L=346)  1150    135        2 ( 410/1637- 512/1735)  76.0 E63-E56 N19-N27
      4  489 1648 (N21E53 L=352)   500     93        5 ( 450/1614- 537/1689)  51.3 E56-E50 N18-N24
      5  765 1582 (N21E29 L= 17)   402    201        0 ( 685/1550- 836/1612)   5.2 E34-E23 N18-N22
      6 1536 1460 (N31W24 L= 69)  2151    474        1 (1345/1298-1664/1574) -18.4 W08-W35 N24-N37
      7 1124 1426 (N18E03 L= 43)   532    230        1 (1053/1345-1198/1507) -28.6 E08-W03 N15-N22
      8  693 1236 (S01E24 L= 22)   372     99        0 ( 660/1208- 718/1285)  76.0 E25-E22 S03-N01
      9 1827 1077 (N16W53 L= 99)  3065    472        9 (1750/ 877-1881/1256) -53.1 W43-W59 N04-N24
     10  523  922 (S23E29 L= 16)   348     90        0 ( 511/ 876- 536/ 967) -80.5 E30-E28 S25-S19
     11  163  823 (S33E68 L=337)  1040    162        1 ( 135/ 745- 209/ 883)  25.3 E79-E59 S38-S29
     12 1720  648 (S11W56 L=102)   382     38        0 (1704/ 619-1736/ 672) -56.3 W54-W58 S12-S09
     13  539  622 (S41E26 L= 20)   272     73        0 ( 502/ 603- 581/ 637)  29.7 E29-E22 S43-S39
     14 1428  510 (S27W40 L= 86)   847    225        2 (1339/ 474-1533/ 571)   0.0 W31-W50 S28-S24
     15 1493  314 (S36W62 L=107)  1310    108        2 (1440/ 275-1546/ 347)   8.7 W53-W68 S38-S33
     16 1391  273 (S41W57 L=103)   460     42        0 (1370/ 245-1410/ 301) -11.3 W52-W62 S42-S39
  

    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