Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  353 1495 (N17E58 L=283)  1895    123        4 ( 261/1377- 431/1638)  78.1 E65-E54 N08-N27
      2  382 1431 (N13E53 L=289)   420     72        0 ( 353/1395- 412/1467)  90.0 E54-E50 N10-N15
      3 1623 1223 (N22W37 L= 18)   272     67        0 (1593/1186-1647/1257) -45.0 W34-W38 N19-N23
      4 1248 1173 (N12W11 L=352)   306     81        0 (1235/1129-1265/1221) -68.2 W09-W12 N09-N14
      5  157 1159 (S06E69 L=272)  1074    204        1 ( 121/1119- 183/1193) -13.0 E78-E64 S08-S04
      6  730 1162 (N03E20 L=321)   209     28        2 ( 705/1142- 753/1184)   0.0 E21-E18 N01-N04
      7 1291 1097 (N08W15 L=356)   258     82        1 (1254/1083-1334/1119)  11.3 W12-W17 N07-N08
      8  155 1021 (S15E66 L=275)   887    120        2 ( 136/ 974- 181/1056)  45.0 E68-E62 S17-S12
      9 1599  536 (S20W50 L= 31)   398    121        0 (1568/ 513-1635/ 564)  -7.1 W46-W54 S21-S18
     10 1245  422 (S34W28 L= 10)   269     65        0 (1213/ 408-1280/ 435)  26.6 W26-W30 S35-S33
  

    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