Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  421 1508 (N27E46 L=124)   782    163        0 ( 406/1432- 442/1592) -54.0 E50-E41 N22-N33
      2 1177 1310 (N14W09 L=178)   281     83        0 (1133/1294-1221/1327) -11.3 W06-W11 N13-N15
      3  634 1198 (N07E24 L=145)   223     49        1 ( 618/1164- 648/1229) -76.0 E25-E23 N04-N08
      4  715 1164 (N05E19 L=150)   246     72        1 ( 687/1128- 749/1185) -14.0 E20-E16 N02-N06
      5 1293  818 (S16W17 L=186)  1167    421        3 (1162/ 668-1373/ 983)  57.7 W09-W21 S25-S05
      6 1719  756 (S18W50 L=219)   831    175        5 (1667/ 703-1794/ 804)  31.0 W45-W56 S21-S15
      7  152  725 (S20E75 L= 95)  2064    129        4 ( 124/ 658- 185/ 785) -76.0 E82-E69 S23-S15
      8 1258  474 (S38W18 L=188)  1106    346        1 (1104/ 455-1432/ 494)   4.4 W06-W32 S39-S36
      9 1089  299 (S52W07 L=176)   569    113        0 (1028/ 282-1142/ 317)  -5.2 W01-W12 S53-S50
  

    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