Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1043 1517 (N36W14 L=296)  1380    313        4 ( 912/1488-1201/1541) -19.3 W04-W24 N31-N39
      2 1427 1483 (N24W38 L=320)  1061    254        4 (1318/1403-1528/1530) -53.1 W32-W42 N16-N28
      3 1850 1184 (S05W62 L=344)  2313    237        7 (1788/1115-1918/1254) -37.6 W56-W72 S10-N00
      4  790 1091 (N16E13 L=269)   703    188        3 ( 762/ 984- 823/1162)  84.8 E14-E09 N09-N20
      5  330 1033 (N20E47 L=234)   682    172        0 ( 294/ 981- 392/1101) -45.0 E50-E42 N16-N24
      6  315  878 (N11E51 L=231)   444    147        1 ( 297/ 821- 337/ 956) -90.0 E52-E49 N07-N16
      7 1760  788 (S26W46 L=327)  1327    269        1 (1676/ 689-1820/ 884)  28.1 W37-W52 S30-S21
      8  412  709 (S00E47 L=235)  1065    212        1 ( 368/ 635- 502/ 756)  12.1 E52-E38 S04-N02
      9  339  532 (S10E61 L=220)  7591    349       14 ( 243/ 389- 439/ 664)  65.8 E74-E52 S21-S01
     10  570  631 (S07E37 L=245)   293     51        2 ( 549/ 600- 590/ 652)  45.0 E38-E34 S09-S05
     11 1283  395 (S39W02 L=284)   270     49        0 (1263/ 369-1304/ 421)  45.0 W00-W03 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