Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1319 1806 (N54W03 L=  1)   549    137        2 (1278/1779-1358/1829)  53.1 W00-W06 N49-N57
      2 1341 1635 (N38W08 L=  6)   272     78        1 (1324/1601-1353/1669) -68.2 W06-W09 N35-N41
      3 1684 1442 (N34W37 L= 35)   323     77        2 (1665/1419-1708/1462)  51.3 W35-W39 N31-N36
      4  635 1422 (N09E31 L=326)   268    130        1 ( 576/1387- 678/1446) -14.0 E35-E27 N07-N10
      5 1655 1187 (N16W36 L= 34)   251     98        2 (1612/1151-1690/1217)  -8.1 W32-W39 N14-N17
      6  989 1137 (S01E04 L=354)   254    135        0 ( 933/1082-1019/1167) -31.0 E07-E02 S04-N00
      7  684  987 (S15E20 L=338)   299     72        0 ( 666/ 949- 701/1031) -78.7 E20-E17 S17-S12
      8  330  775 (S33E44 L=313)  2090    574        3 ( 180/ 683- 481/ 916)  32.7 E60-E32 S40-S22
      9  420  575 (S46E38 L=320)   401    141        0 ( 378/ 539- 451/ 612)  40.6 E43-E34 S48-S42
     10 1362  245 (S45W52 L= 49)   498    106        2 (1331/ 232-1399/ 268)   0.0 W47-W54 S46-S43
     11  861  200 (S67W17 L= 14)   465     80        1 ( 838/ 181- 889/ 219)  68.2 W13-W19 S69-S64
  

    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