Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1237 1403 (N27W18 L=203)   233     71        0 (1225/1367-1249/1440)  68.2 W16-W18 N24-N29
      2 1120 1374 (N26W09 L=195)   445    142        2 (1070/1338-1166/1404) -33.7 W05-W11 N23-N28
      3  374 1070 (N11E44 L=141)  1161    244        2 ( 349/ 959- 404/1185) -70.3 E47-E41 N03-N17
      4 1200  969 (N00W10 L=196)   209     59        1 (1176/ 948-1231/ 994)  14.0 W08-W12 N00-N01
      5 1675  944 (S06W43 L=229)  1005    153        2 (1612/ 903-1755/ 974)  -9.5 W38-W50 S07-S03
      6 1855  899 (S11W63 L=248)   563     79        1 (1834/ 862-1876/ 937) -59.0 W59-W66 S13-S08
      7  498  753 (S09E37 L=149)   321     83        1 ( 464/ 730- 531/ 786) -36.9 E38-E34 S10-S07
      8  357  694 (S13E50 L=135)   355     86        1 ( 329/ 662- 391/ 730) -90.0 E52-E46 S14-S10
      9  992  548 (S25E06 L=180)  1669    495        4 ( 861/ 377-1147/ 700)  43.7 E17-W05 S37-S15
     10 1766  603 (S30W59 L=245)   527     64        1 (1732/ 580-1794/ 629)  14.0 W54-W62 S31-S28
     11 1491  581 (S27W30 L=216)   547    130        0 (1447/ 532-1545/ 626)  35.5 W26-W34 S30-S24
     12 1645  558 (S31W45 L=231)   673    151        0 (1606/ 493-1673/ 624)  77.5 W43-W46 S35-S26
  

    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