Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  553 1778 (N26E59 L= 90)   617    144        1 ( 534/1737- 571/1814) -56.3 E62-E56 N23-N29
      2  873 1689 (N28E31 L=119)   371     91        0 ( 831/1660- 911/1707)   0.0 E33-E27 N27-N29
      3  492 1592 (N13E51 L= 98)   472     93        1 ( 463/1550- 530/1631)  80.5 E53-E49 N09-N15
      4  993 1504 (N20E16 L=134)   300     96        1 ( 972/1459-1022/1548) -45.0 E17-E12 N17-N22
      5 1678 1092 (N17W37 L=187)   523    122        0 (1636/1059-1720/1118)  -8.1 W34-W41 N14-N18
      6 1138  636 (S24W18 L=168)   217     75        0 (1094/ 622-1181/ 647)  18.4 W15-W21 S25-S23
      7 1065  484 (S35W21 L=170)   252     56        0 (1044/ 460-1093/ 511) -11.3 W18-W23 S36-S34
      8 1273  446 (S31W35 L=185)   270     76        1 (1253/ 423-1300/ 479)  -9.5 W32-W38 S31-S29
      9 1636  406 (S20W64 L=213)   719    189        2 (1606/ 373-1666/ 440)   0.0 W59-W69 S20-S18
  

    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