Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  802 1491 (N39E09 L=127)   594     87        2 ( 757/1463- 844/1521)   8.1 E12-E05 N36-N40
      2  249 1262 (N29E59 L= 77)   516     57        1 ( 234/1232- 264/1296) -63.4 E62-E56 N26-N30
      3  947 1199 (N18E02 L=134)   613    220        0 ( 877/1115- 991/1266) -68.2 E06-W00 N12-N22
      4  423 1198 (N24E40 L= 96)   270     70        0 ( 412/1158- 432/1241) -78.7 E40-E38 N21-N26
      5  852 1180 (N18E09 L=128)   324    119        1 ( 833/1119- 879/1234) -90.0 E09-E07 N13-N21
      6  247 1179 (N23E58 L= 78)   569     95        5 ( 210/1141- 277/1211) -18.4 E63-E53 N21-N25
      7 1879  950 (S13W64 L=201)  1237    158        1 (1851/ 880-1908/1014)  37.9 W60-W69 S17-S10
      8  701  949 (N06E21 L=115)   399    150        0 ( 663/ 887- 744/1013)  45.0 E24-E17 N02-N09
      9  624  535 (S20E34 L=102)   509    115        1 ( 588/ 492- 673/ 580)  29.1 E38-E29 S22-S16
     10  407  461 (S23E58 L= 79)   995    154        1 ( 351/ 428- 470/ 495)   3.2 E69-E51 S25-S20
     11  508  381 (S30E53 L= 84)  1031    138        0 ( 449/ 331- 565/ 427)  15.3 E66-E44 S33-S26
     12  922  341 (S38E20 L=117)   575    125        1 ( 849/ 320- 979/ 368)   0.0 E25-E14 S39-S36
  

    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