Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  464 1597 (N38E49 L=238)   482    116        1 ( 437/1573- 495/1623) -15.9 E52-E45 N36-N40
      2  685 1202 (N13E22 L=264)   922    349        1 ( 562/1142- 785/1303)  31.6 E30-E16 N09-N19
      3  229 1268 (N15E62 L=225)   425     69        1 ( 204/1250- 251/1284)   0.0 E64-E58 N13-N15
      4  625  879 (S07E25 L=262)   367    134        0 ( 577/ 851- 676/ 904) -15.9 E28-E21 S08-S05
      5 1841  813 (S10W63 L=350)   785    161        0 (1790/ 775-1874/ 849)  26.6 W57-W67 S12-S07
      6  927  632 (S22E06 L=281)  1238    383        4 ( 807/ 553-1016/ 775)  49.4 E13-W00 S27-S12
      7  664  705 (S18E23 L=264)   222     35        1 ( 640/ 684- 685/ 729) -56.3 E24-E21 S19-S16
      8 1713  585 (S25W54 L=341)   382     50        0 (1688/ 564-1730/ 611)  33.7 W51-W55 S26-S22
      9  413  400 (S42E53 L=234)  1400    291        1 ( 364/ 373- 466/ 436)  15.9 E61-E47 S45-S38
  

    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