Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1053 1627 (N34E01 L=259)   281     53        1 (1029/1602-1084/1645)  56.3 E03-W00 N31-N35
      2 1650 1577 (N35W47 L=307)   517     92        0 (1607/1555-1689/1595)  11.3 W41-W51 N33-N36
      3 1546 1547 (N32W35 L=296)   485     89        1 (1531/1498-1564/1591)  66.8 W33-W37 N28-N35
      4  285 1492 (N22E60 L=200)  1431    121        2 ( 250/1437- 323/1546)  74.1 E64-E54 N18-N26
      5 1831 1454 (N29W65 L=325)   969    101        1 (1790/1423-1877/1476)  -9.5 W58-W72 N27-N30
      6  250 1362 (N14E58 L=202)   843    108        1 ( 195/1334- 290/1402)  16.7 E64-E54 N12-N16
      7  433 1338 (N12E40 L=220)   278    142        0 ( 392/1290- 495/1364)  39.8 E43-E36 N09-N14
      8  630 1275 (N09E26 L=235)   395    145        0 ( 573/1255- 703/1299)  -6.3 E29-E20 N07-N10
      9  163  920 (S12E64 L=197)  1677    267        6 ( 126/ 782- 197/1033) -90.0 E68-E58 S20-S05
     10   89  970 (S09E76 L=185)   916     73        2 (  68/ 931- 113/ 999)  18.4 E82-E70 S10-S07
     11 1520  598 (S27W38 L=298)  1632    271        3 (1430/ 492-1635/ 674)  55.3 W32-W45 S34-S21
     12 1321  576 (S30W23 L=284)   437    137        0 (1260/ 542-1375/ 618)  40.6 W19-W26 S32-S26
     13  642  578 (S33E25 L=235)   865    296        0 ( 513/ 504- 742/ 614) -13.2 E35-E18 S38-S30
     14 1647  523 (S30W52 L=313)   410     96        0 (1606/ 495-1681/ 560)  51.3 W49-W54 S32-S27
     15  457  518 (S38E43 L=217)   507    105        3 ( 414/ 487- 526/ 546)  12.1 E49-E35 S40-S35
     16 1604  328 (S44W65 L=326)  1008    100        1 (1561/ 299-1652/ 360)  -4.4 W56-W82 S44-S41
     17 1012  318 (S51W04 L=265)  1387    294        2 ( 858/ 294-1132/ 355)   6.8 E10-W15 S52-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