Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  924 1510 (N18E20 L= 82)   719    204        0 ( 897/1420- 949/1605) -63.4 E22-E16 N12-N24
      2  321 1567 (N08E68 L= 34)  1302    141        4 ( 285/1525- 360/1600)  39.8 E72-E63 N05-N10
      3 1297 1509 (N29W02 L=104)   614    103        2 (1262/1447-1323/1559) -50.2 E01-W05 N25-N31
      4 1113 1374 (N15E05 L= 97)   566    178        0 (1078/1305-1130/1464) -48.8 E09-E02 N10-N19
      5  520 1367 (S00E39 L= 62)   422    111        2 ( 482/1323- 560/1408) -15.9 E43-E35 S02-N01
      6 1643 1142 (N18W34 L=135)   270     82        1 (1617/1109-1669/1171)  78.7 W31-W34 N15-N20
      7  899 1097 (S06E09 L= 93)   501     93        1 ( 865/1043- 925/1154) -81.9 E09-E07 S09-S02
      8  311  911 (S31E45 L= 57)  3135    634        7 ( 158/ 760- 437/1098)  50.0 E60-E36 S42-S17
      9 1652  967 (N09W39 L=140)   333     89        1 (1620/ 944-1690/ 986)  31.0 W36-W41 N06-N10
     10 1867  937 (N15W59 L=161)   483     66        2 (1844/ 906-1886/ 966) -71.6 W56-W60 N13-N16
     11 1696  829 (N03W47 L=148)  1162    381        1 (1561/ 790-1814/ 881)  24.2 W35-W55 S01-N09
     12  640  769 (S31E16 L= 85)   809    257        0 ( 567/ 672- 712/ 867)  68.2 E20-E13 S38-S23
     13  974  784 (S21W04 L=106)   752    205        1 ( 871/ 719-1057/ 841)  57.5 W00-W07 S27-S15
     14 1199  707 (S19W19 L=121)   551    152        0 (1128/ 666-1286/ 737)  48.8 W16-W23 S22-S14
     15 1693  709 (S03W51 L=153)   339     57        1 (1668/ 683-1718/ 733)  63.4 W49-W53 S05-S01
     16 1564  556 (S15W48 L=150)   400     78        0 (1527/ 535-1597/ 582)  53.1 W45-W50 S16-S12
     17 1295  435 (S31W37 L=139)  1339    249        3 (1187/ 393-1416/ 474)  37.6 W31-W44 S35-S25
     18 1548  412 (S23W56 L=158)   468     74        1 (1512/ 382-1580/ 443)  68.2 W54-W58 S25-S20
  

    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