Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1103 1740 (N53W22 L=229)   753    111        2 (1066/1692-1128/1791)  33.7 W15-W27 N49-N57
      2  338 1282 (N28E50 L=158)   431    114        2 ( 310/1226- 353/1327) -80.5 E53-E47 N24-N30
      3  146 1162 (N20E74 L=133)   799     72        2 ( 116/1139- 174/1183)  13.0 E81-E68 N18-N21
      4  906 1114 (N13E06 L=201)   227     86        1 ( 875/1079- 934/1145) -63.4 E07-E04 N10-N14
      5  997 1021 (N06E02 L=206)   252    100        2 ( 952/1002-1043/1042)   0.0 E04-W01 N05-N06
      6  387  903 (N04E44 L=163)   276     53        1 ( 368/ 883- 410/ 935)  76.0 E46-E42 N02-N06
      7  569  478 (S25E40 L=167)   311     60        2 ( 538/ 452- 591/ 501)  18.4 E43-E37 S26-S23
      8  735  409 (S31E30 L=178)   656    159        3 ( 672/ 371- 808/ 442)  13.0 E36-E23 S34-S28
      9 1404  421 (S40W22 L=229)   306     73        0 (1364/ 411-1448/ 431)  -8.1 W18-W25 S40-S38
     10 1207  310 (S46W04 L=212)   996    185        2 (1132/ 277-1280/ 370)  17.1 E02-W11 S49-S41
  

    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