Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1271 1665 (N39W05 L=237)  1305    360        1 (1120/1589-1421/1749) -13.5 E08-W17 N36-N43
      2  386 1635 (N23E62 L=170)  1197     95        4 ( 352/1586- 408/1678) -49.4 E67-E57 N19-N26
      3  511 1121 (S08E33 L=199)   228     61        0 ( 493/1088- 536/1151)  63.4 E34-E31 S09-S05
      4  525  903 (S21E30 L=202)   524    183        0 ( 484/ 838- 552/ 971)  90.0 E32-E28 S25-S16
      5 1228  866 (S12W15 L=247)   214     96        0 (1193/ 840-1279/ 894)  51.3 W13-W17 S14-S09
      6  147  834 (S28E70 L=162)   813     47        0 ( 112/ 818- 177/ 852)   0.0 E79-E65 S28-S26
      7  751  774 (S25E13 L=219)   190     99        0 ( 708/ 746- 781/ 811)  51.3 E15-E11 S27-S22
      8 1270  655 (S24W22 L=254)  2890    824        3 ( 960/ 544-1654/ 802)  35.0 W04-W44 S35-S07
      9 1278  772 (S17W20 L=252)   261     78        0 (1233/ 754-1314/ 791)  36.9 W17-W21 S18-S15
     10 1661  628 (S17W52 L=284)  1241    225        4 (1546/ 591-1739/ 675)  21.8 W42-W57 S20-S12
     11  840  580 (S36E05 L=228)   405    115        1 ( 793/ 549- 891/ 615)  -7.1 E08-W00 S38-S34
  

    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