Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  591 1663 (N21E48 L=115)  1681    342        2 ( 559/1542- 635/1778) -57.1 E54-E42 N12-N29
      2  370 1533 (N08E60 L=103)  1740    459        4 ( 344/1380- 397/1684) -48.4 E70-E51 N00-N18
      3  236 1378 (S03E67 L= 97)   751    100        1 ( 210/1340- 269/1414)  63.4 E70-E63 S06-N00
      4  735 1380 (N07E26 L=137)   375     84        0 ( 687/1352- 775/1412)  63.4 E28-E24 N04-N10
      5 1116 1313 (N14E03 L=161)   299    134        0 (1068/1249-1157/1375) -32.0 E07-W01 N10-N15
      6  887 1343 (N09E16 L=147)   234     86        0 ( 841/1328- 925/1368)  45.0 E18-E14 N07-N11
      7 1742 1149 (N23W42 L=206)   543     86        1 (1723/1102-1756/1198) -59.0 W40-W44 N20-N26
      8  359 1044 (S21E43 L=121)   535    148        2 ( 319/ 999- 416/1076)  51.3 E46-E38 S23-S18
      9  185 1041 (S24E59 L=104)   418     89        1 ( 163/1008- 210/1071)  53.1 E61-E56 S26-S22
     10  618  863 (S25E20 L=143)  1285    459        1 ( 445/ 811- 834/ 949)  36.9 E32-E08 S32-S14
     11 1133  783 (S15W13 L=176)   375    162        1 (1051/ 754-1202/ 804)  49.4 W09-W16 S19-S12
     12 1839  737 (N04W66 L=229)   782     93        1 (1821/ 688-1855/ 778) -51.3 W62-W69 N00-N06
     13 1169  684 (S20W18 L=182)   861    265        1 (1056/ 660-1315/ 711)  21.8 W11-W26 S23-S14
     14 1518  549 (S17W46 L=209)   506    114        3 (1475/ 519-1559/ 578)  51.3 W42-W48 S19-S14
     15  994  492 (S36W16 L=179)   318     82        0 ( 958/ 472-1029/ 523)  68.2 W13-W17 S38-S33
     16 1581  324 (S26W66 L=229)   653    136        2 (1546/ 307-1620/ 339)  23.2 W62-W69 S27-S24
     17 1457  291 (S32W58 L=221)   532     91        3 (1432/ 268-1484/ 314)   9.5 W54-W60 S34-S31
  

    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