Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  688 1834 (N34E57 L=175)  1063    116        4 ( 649/1786- 722/1877) -17.5 E69-E50 N30-N37
      2  518 1696 (N20E58 L=175)   472     73        1 ( 486/1677- 546/1713)  45.0 E60-E53 N18-N21
      3  448 1616 (N13E58 L=175)  2486    296        7 ( 337/1559- 548/1664)  17.5 E71-E49 N08-N17
      4  616 1626 (N17E44 L=188)   545    118        1 ( 573/1575- 662/1663)  63.4 E47-E41 N13-N19
      5 1284 1521 (N29W01 L=233)   227     96        0 (1262/1477-1297/1561) -59.0 E01-W02 N26-N31
      6  676 1500 (N11E34 L=198)   433    125        0 ( 639/1457- 701/1558) -74.1 E36-E31 N07-N14
      7  556 1523 (N10E43 L=189)   273     46        1 ( 536/1501- 580/1543) -14.0 E45-E41 N08-N10
      8 1700 1439 (N40W35 L=268)   567    135        0 (1657/1380-1732/1490) -35.5 W30-W38 N36-N42
      9 1647 1323 (N29W31 L=263)   339     92        0 (1633/1278-1667/1367) -51.3 W29-W33 N26-N31
     10  382 1162 (S15E44 L=188)   424     95        4 ( 353/1131- 411/1191)  76.0 E46-E41 S16-S12
     11  486 1135 (S14E35 L=197)   250     62        1 ( 472/1106- 509/1158)  71.6 E36-E33 S15-S12
     12 1911  940 (N17W67 L=299)   672     87        0 (1874/ 916-1939/ 973)  33.7 W61-W71 N14-N20
     13  528  892 (S27E27 L=205)   340    115        1 ( 472/ 870- 579/ 913)   7.1 E31-E23 S27-S26
     14  286  838 (S36E49 L=183)   332     98        0 ( 248/ 809- 328/ 861)  26.6 E54-E44 S38-S33
     15 1652  500 (S15W60 L=293)  1811    374        1 (1562/ 390-1742/ 638) -83.7 W55-W65 S23-S05
     16 1142  570 (S29W21 L=254)   241    103        0 (1103/ 554-1198/ 591)   8.1 W18-W25 S29-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