Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  655 1678 (N40E41 L=124)   611    109        1 ( 616/1635- 690/1718) -31.0 E48-E37 N36-N42
      2  436 1611 (N32E59 L=106)   550     73        0 ( 418/1574- 457/1650) -24.4 E64-E53 N29-N34
      3 1287 1582 (N41W14 L=179)   760    152        2 (1227/1543-1356/1624) -15.3 W08-W19 N38-N43
      4  652 1427 (N22E30 L=135)   586    132        1 ( 622/1372- 685/1486) -76.0 E31-E28 N18-N26
      5 1511 1214 (N18W30 L=195)   310    114        0 (1486/1167-1549/1266) -45.0 W28-W33 N15-N20
      6  307  890 (S15E50 L=115)  4132    550        8 ( 193/ 658- 426/1082)  41.4 E74-E40 S31-S01
      7 1705  730 (S10W52 L=217)   586    110        0 (1690/ 680-1719/ 786) -81.9 W50-W53 S13-S06
      8 1666  716 (S12W48 L=214)   336    129        0 (1627/ 675-1686/ 771)  80.5 W46-W50 S14-S08
      9 1258  744 (S14W18 L=183)   317     87        0 (1213/ 728-1301/ 761)   0.0 W14-W20 S14-S12
     10 1240  556 (S26W20 L=185)   372     88        0 (1193/ 524-1281/ 584)  38.7 W17-W22 S28-S23
     11 1632  504 (S26W55 L=220)  1461    154        3 (1556/ 465-1703/ 546)  21.0 W48-W62 S29-S22
     12 1278  435 (S34W27 L=192)   655    115        1 (1225/ 416-1334/ 458)  14.0 W22-W30 S35-S32
  

    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