Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1305 1507 (N32W03 L=123)   370    117        0 (1273/1452-1333/1549) -33.7 W00-W06 N28-N33
      2 1722 1400 (N40W41 L=160)   452     78        0 (1690/1378-1756/1420)   0.0 W36-W44 N37-N40
      3 1641 1272 (N28W33 L=153)  1298    242        1 (1572/1164-1712/1334) -23.2 W26-W40 N22-N30
      4  562 1289 (S01E34 L= 85)   690    134        0 ( 495/1241- 606/1322)  -6.3 E39-E30 S03-N00
      5 1468 1221 (N19W21 L=141)   721    271        0 (1365/1139-1537/1321)  66.4 W16-W23 N11-N27
      6 1788 1243 (N31W48 L=168)   439     88        0 (1776/1189-1800/1288)  81.9 W46-W49 N27-N34
      7  335 1132 (S15E48 L= 72)   321    154        0 ( 290/1094- 406/1152)  33.7 E51-E41 S18-S12
      8 1811  993 (N17W54 L=174)   805    123        1 (1782/ 954-1837/1053) -33.7 W49-W56 N14-N19
      9  179  914 (S32E64 L= 56)  2767    245        4 ( 127/ 815- 223/1019)  33.7 E80-E57 S38-S24
     10 1356  791 (S08W25 L=145)  1492    470        1 (1159/ 727-1518/ 850)  36.4 W15-W34 S15-S01
     11 1662  782 (N00W46 L=166)   671    143        2 (1637/ 709-1684/ 850) -76.0 W44-W48 S04-N04
     12  192  801 (S39E64 L= 56)   491     47        1 ( 167/ 782- 214/ 824)  16.7 E69-E59 S40-S36
     13  731  708 (S30E09 L=111)   676    237        1 ( 656/ 624- 824/ 768)  56.3 E12-E04 S37-S24
     14 1717  548 (S10W64 L=184)   509     18        0 (1694/ 528-1741/ 566)   9.5 W60-W67 S11-S08
     15 1352  313 (S35W52 L=172)   633    128        1 (1308/ 299-1420/ 335)  16.7 W46-W56 S37-S32
     16 1531  299 (S29W70 L=190)   999     85        2 (1497/ 273-1568/ 332)   3.6 W62-W80 S30-S27
  

    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