Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  611 1745 (N30E52 L=355)  1396    281        4 ( 555/1676- 649/1801) -29.7 E59-E44 N25-N33
      2 1058 1538 (N28E12 L= 35)  1331    364        2 ( 910/1418-1204/1662) -15.1 E26-W00 N22-N32
      3  422 1580 (N15E57 L=350)  1935    393        4 ( 341/1467- 550/1661)  57.3 E63-E48 N06-N21
      4  538 1499 (N12E43 L=  4)  1116    368        0 ( 495/1333- 575/1624) -69.8 E48-E39 N01-N20
      5 1479 1498 (N37W18 L= 65)   559    123        1 (1419/1457-1529/1536) -11.3 W12-W22 N35-N39
      6  669 1430 (N11E32 L= 15)   262     78        0 ( 652/1389- 684/1464) -36.9 E34-E29 N09-N13
      7 1638 1348 (N32W33 L= 80)   315     80        0 (1600/1324-1671/1365)   0.0 W29-W35 N31-N33
      8  575 1238 (S02E32 L= 15)   470    116        1 ( 528/1220- 633/1267)   0.0 E35-E27 S03-N00
      9  721  564 (S39E08 L= 40)   422     91        0 ( 677/ 537- 772/ 590)   0.0 E10-E03 S40-S37
  

    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