Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  634 1331 (N30E22 L=207)   221     89        1 ( 592/1317- 674/1344) -18.4 E25-E19 N29-N31
      2 1220 1217 (N15W15 L=245)   541    175        2 (1135/1187-1297/1259) -33.7 W09-W19 N11-N18
      3  173 1195 (N25E64 L=166)  1129    299        1 ( 117/1171- 230/1220)  13.0 E70-E57 N23-N26
      4  458 1072 (N16E37 L=193)  1439    311        2 ( 413/ 933- 523/1206) -70.6 E39-E32 N07-N24
      5  557 1171 (N21E29 L=201)   271     61        0 ( 530/1144- 588/1201) -45.0 E30-E26 N19-N23
      6 1181  967 (N01W08 L=238)   289     89        1 (1128/ 935-1219/ 998) -38.7 W05-W10 S01-N03
      7  237  896 (N07E59 L=171)   436     96        2 ( 211/ 871- 269/ 925) -18.4 E61-E55 N05-N08
      8  279  843 (N03E55 L=175)   388     73        2 ( 264/ 810- 293/ 876)  53.1 E56-E53 N01-N05
      9 1040  542 (S24E06 L=223)   308    113        0 (1024/ 487-1053/ 589)  74.1 E07-E04 S27-S20
     10 1433  261 (S58W26 L=255)   486     86        0 (1395/ 245-1474/ 274)  -5.7 W20-W31 S59-S56
  

    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