Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1076 1902 (N58W57 L=314)   820     78        1 (1039/1880-1105/1926) -36.9 W51-W66 N55-N60
      2 1610  891 (S18W32 L=289)   260     62        0 (1595/ 853-1624/ 929) -90.0 W30-W32 S20-S15
      3 1702  839 (S25W38 L=295)   695    143        0 (1636/ 793-1764/ 879)   0.0 W32-W44 S25-S23
      4 1448  771 (S20W18 L=275)   223     89        1 (1416/ 736-1485/ 799)   0.0 W15-W20 S21-S18
      5 1524  637 (S31W21 L=278)   820    206        1 (1471/ 604-1595/ 713) -55.0 W16-W26 S35-S25
      6 1427  632 (S28W14 L=270)   432    120        0 (1386/ 580-1450/ 680)  39.8 W09-W16 S30-S25
      7 1226  563 (S27E01 L=256)   438    104        0 (1200/ 508-1255/ 603)  45.0 E04-W01 S29-S24
      8 1768  522 (S51W51 L=308)   746     69        0 (1733/ 509-1799/ 539) -11.3 W43-W61 S54-S48
  

    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