Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1295 1703 (N44E06 L=196)   547    153        0 (1236/1642-1339/1757) -26.6 E12-W00 N40-N46
      2  444 1398 (S00E47 L=155)   775    108        1 ( 415/1331- 467/1443) -26.6 E50-E42 S03-N02
      3 1365 1283 (N18W12 L=214)   738    183        0 (1291/1212-1419/1331) -16.7 W06-W16 N15-N19
      4  944  830 (S19W01 L=203)  1247    464        0 ( 736/ 744-1139/ 874)  43.5 E09-W10 S29-S11
      5 1541  850 (S00W35 L=237)   233     66        0 (1511/ 831-1580/ 867)   0.0 W32-W38 S01-N00
      6  516  708 (S39E25 L=177)  1123    302        1 ( 473/ 579- 586/ 807)  90.0 E28-E20 S46-S30
      7 1045  670 (S25W12 L=214)   589    174        0 ( 974/ 635-1129/ 700)  56.3 W08-W16 S29-S20
      8 1100  549 (S30W20 L=222)   698    228        0 ( 979/ 533-1194/ 566)  28.3 W12-W25 S34-S26
      9  446  425 (S63E35 L=168)  1319    159        4 ( 381/ 391- 513/ 475)  -5.4 E45-E23 S65-S59
     10 1392  426 (S27W46 L=248)   691    176        0 (1326/ 399-1470/ 463)   7.6 W38-W53 S28-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