Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  600 1572 (N24E42 L=168)  3070    526        5 ( 455/1454- 825/1726) -11.5 E70-E21 N19-N31
      2  974 1436 (N22E12 L=198)   270     86        0 ( 939/1419-1022/1460)   0.0 E14-E08 N21-N23
      3  833 1264 (N09E16 L=193)   832    331        4 ( 718/1217- 945/1342)  -3.8 E24-E09 N06-N11
      4  665 1188 (N02E25 L=185)   391    102        0 ( 618/1166- 719/1207)  18.4 E27-E21 N00-N02
      5 1119 1108 (N05W04 L=213)   449    190        1 (1093/1019-1139/1201) -84.8 W02-W04 N00-N11
      6  521  894 (S18E30 L=180)   436     70        0 ( 499/ 857- 547/ 936)  78.7 E31-E27 S20-S15
      7  572  774 (S25E25 L=185)   632    192        3 ( 501/ 731- 642/ 838) -15.3 E30-E19 S26-S22
      8 1661  469 (S22W63 L=273)   449     54        0 (1631/ 448-1690/ 490)  36.9 W60-W66 S23-S19
      9 1287  368 (S37W37 L=246)   324     62        1 (1265/ 342-1312/ 398) -23.2 W33-W40 S38-S34
     10 1428  372 (S33W48 L=258)   573     85        2 (1382/ 353-1468/ 395)  14.0 W44-W52 S34-S30
  

    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