Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  498 1373 (N11E41 L= 87)   257     55        0 ( 479/1344- 520/1399)  76.0 E42-E39 N09-N13
      2 1047 1347 (N19E04 L=124)   281     53        1 (1022/1329-1081/1369) -14.0 E05-E01 N17-N20
      3 1637 1255 (N24W38 L=166)   410     91        2 (1585/1234-1677/1273)  14.0 W33-W41 N22-N25
      4 1229 1057 (N05W12 L=140)   919    248        1 (1145/ 984-1330/1129) -19.7 W05-W19 N02-N07
      5 1877  957 (N10W66 L=194)   554     84        2 (1850/ 924-1900/ 991) -23.2 W61-W69 N08-N11
      6 1832  868 (N04W62 L=190)   680     94        5 (1787/ 844-1870/ 900)  16.7 W56-W66 N02-N06
      7  851  732 (S21E06 L=122)   497    203        2 ( 824/ 644- 896/ 819) -59.0 E08-E01 S26-S15
      8  620  629 (S32E21 L=107)   265     62        0 ( 605/ 596- 635/ 662) -78.7 E22-E19 S34-S29
      9 1672  549 (S18W58 L=186)   775    126        1 (1637/ 501-1706/ 601)  90.0 W56-W59 S21-S14
     10 1155  536 (S28W18 L=146)   369     90        1 (1125/ 488-1178/ 579)  90.0 W16-W19 S31-S25
     11  855  427 (S42E00 L=128)   823    167        2 ( 787/ 405- 927/ 459)  24.0 E06-W05 S44-S38
  

    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