Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  665 1403 (N33E14 L=178)  1085    283        4 ( 533/1337- 759/1449) -20.6 E24-E08 N28-N37
      2 1551 1331 (N08W39 L=232)   458    110        3 (1520/1295-1581/1363)  76.0 W36-W42 N05-N09
      3  162  968 (N17E62 L=130)  2260    362        6 ( 121/ 796- 198/1116)  78.1 E68-E57 N07-N26
      4  856  747 (S09E16 L=176)  3568    860        9 ( 684/ 490-1005/ 922)  11.8 E29-E06 S22-N03
      5 1284  879 (S11W12 L=204)   321     96        1 (1238/ 862-1335/ 897) -38.7 W09-W14 S12-S08
      6 1323  771 (S18W12 L=204)   295     79        3 (1278/ 743-1366/ 795)  -9.5 W08-W14 S20-S16
      7  681  638 (S11E29 L=163)   318    102        0 ( 646/ 601- 710/ 685)  26.6 E32-E26 S13-S09
      8 1074  566 (S25E08 L=184)   517    150        0 (1011/ 529-1152/ 614)  11.3 E12-E02 S25-S23
      9  517  310 (S29E57 L=135)   692    110        2 ( 493/ 274- 541/ 349)  71.6 E59-E54 S31-S25
  

    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