Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  220 1208 (N32E53 L=179)  1343    208        3 ( 164/1155- 272/1249) -35.0 E58-E47 N27-N35
      2  950 1097 (N09E03 L=229)   180     65        1 ( 910/1076- 975/1121)  11.3 E05-W00 N08-N10
      3  308  660 (S02E57 L=175)  2620    458        7 ( 224/ 570- 446/ 736) -19.3 E69-E47 S08-N02
      4  296  585 (S06E62 L=170)   440     67        0 ( 267/ 559- 319/ 603) -26.6 E64-E59 S08-S05
      5  269  525 (S10E69 L=162)   890    160        2 ( 243/ 474- 301/ 560)  45.0 E72-E63 S13-S07
      6  439  464 (S16E54 L=178)  1361    328        0 ( 401/ 376- 468/ 553)  32.0 E63-E47 S21-S10
      7  578  384 (S24E47 L=184)   384     90        1 ( 557/ 354- 610/ 405)  18.4 E50-E44 S25-S22
  

    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