Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1262 1713 (N44W39 L= 83)   507    125        1 (1232/1680-1311/1737)   6.3 W34-W43 N42-N45
      2 1104 1355 (N25W12 L= 56)   342    107        1 (1065/1320-1156/1379)   8.1 W08-W15 N23-N26
      3  249 1340 (N39E55 L=349)   442     78        1 ( 222/1314- 271/1363)  36.9 E57-E52 N37-N40
      4  746 1330 (N30E12 L= 32)   313    123        1 ( 682/1305- 802/1346) -26.6 E16-E08 N27-N32
      5  287 1086 (N22E50 L=354)   866    206        3 ( 216/1059- 347/1119)  -8.7 E57-E44 N20-N24
      6  525  949 (N11E33 L= 12)   322     68        0 ( 491/ 928- 560/ 967) -26.6 E34-E29 N09-N12
      7  141  930 (N13E70 L=334)   873    101        0 ( 116/ 887- 163/ 965) -59.0 E73-E67 N10-N15
      8  391  453 (S18E59 L=345)   807    174        0 ( 359/ 387- 422/ 523)  78.7 E61-E56 S23-S13
      9 1329  344 (S45W07 L= 51)   392     74        0 (1303/ 328-1365/ 365) -45.0 W04-W09 S47-S42
  

    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