Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1034 1629 (N39W20 L=284)   555    139        0 ( 971/1607-1111/1645) -15.3 W14-W25 N37-N41
      2 1633 1569 (N18W59 L=324)  2072    355        2 (1484/1465-1776/1633) -46.7 W50-W67 N08-N25
      3 1407 1496 (N21W37 L=302)   887    236        0 (1288/1456-1505/1544) -42.3 W31-W42 N16-N26
      4 1232 1470 (N24W25 L=290)   367    128        0 (1205/1413-1272/1527)  82.9 W23-W26 N19-N27
      5 1367 1352 (N14W29 L=293)   280    118        0 (1319/1321-1420/1373) -56.3 W26-W31 N10-N16
      6  325 1254 (N34E44 L=220)  2244    508        4 ( 127/1197- 553/1284) -14.5 E78-E24 N24-N39
      7 1770  803 (S29W46 L=311)  2035    277        1 (1679/ 686-1864/ 892)  21.8 W36-W57 S33-S25
      8 1628  774 (S26W32 L=297)   358    128        0 (1599/ 718-1658/ 829)  45.0 W29-W35 S28-S23
      9  240  661 (N00E65 L=199)   476     56        0 ( 207/ 648- 275/ 674)   0.0 E70-E60 S01-N01
     10 1458  591 (S33W17 L=281)   500    151        0 (1405/ 551-1512/ 661) -59.0 W13-W20 S36-S26
     11  421  596 (S07E50 L=215)   484    132        0 ( 384/ 536- 448/ 648) -90.0 E51-E48 S10-S03
     12  358  560 (S08E58 L=207)   470     94        2 ( 327/ 515- 379/ 597)  78.7 E60-E55 S10-S05
     13  317  536 (S09E64 L=201)   687     78        4 ( 282/ 500- 348/ 564) -53.1 E67-E60 S10-S06
     14  673  467 (S20E37 L=227)   315     77        0 ( 642/ 433- 696/ 509)  33.7 E40-E34 S21-S17
     15 1190  358 (S41E09 L=256)   351     76        2 (1162/ 334-1215/ 397)  18.4 E11-E05 S43-S38
     16  872  141 (S50E53 L=212)   708     93        1 ( 831/ 116- 917/ 161)   3.2 E64-E46 S51-S47
  

    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