Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1199 1718 (N42E13 L=345)   570    107        1 (1168/1672-1226/1755) -26.6 E17-E09 N38-N44
      2  734 1537 (N15E33 L=326)   448    223        0 ( 718/1439- 753/1634) -56.3 E36-E28 N09-N21
      3  850 1557 (N19E26 L=332)   589    190        1 ( 816/1465- 871/1634) -42.0 E32-E22 N14-N23
      4 1650 1386 (N35W31 L= 29)   847    141        2 (1621/1313-1675/1447)  90.0 W27-W33 N29-N39
      5  974 1324 (N09E11 L=347)   886    261        1 ( 895/1248-1076/1396)  65.2 E14-E07 N02-N15
      6 1434 1167 (N14W19 L= 17)   315    184        0 (1376/1104-1500/1233) -15.3 W13-W24 N12-N15
      7  313 1203 (S14E51 L=307)   645    173        1 ( 247/1167- 388/1230)   8.1 E57-E43 S16-S12
      8 1160 1197 (N07W03 L=  1)   219     22        0 (1141/1175-1180/1221) -45.0 W01-W04 N05-N08
      9 1183 1073 (N01W07 L=  6)   774    218        1 (1109/1041-1265/1132)  24.0 W02-W11 S01-N05
     10  774 1009 (S14E14 L=345)   316    118        0 ( 718/ 973- 822/1040)  54.5 E16-E11 S17-S10
     11 1651  940 (N08W40 L= 38)   314     96        1 (1638/ 887-1662/ 988) -51.3 W37-W41 N05-N10
     12  649  963 (S19E20 L=338)   204     44        0 ( 618/ 947- 675/ 977)   0.0 E22-E18 S20-S18
     13  292  822 (S37E48 L=310)  1132    270        1 ( 187/ 769- 414/ 873)  21.8 E65-E35 S42-S30
     14 1134  650 (S23W18 L= 16)  1664    517        4 ( 934/ 573-1363/ 761)  52.0 W08-W27 S33-S10
     15 1803  693 (N01W64 L= 62)   469     97        2 (1776/ 654-1822/ 734)  -7.1 W58-W66 S01-N02
     16  400  537 (S55E38 L=321)   831    120        2 ( 350/ 511- 450/ 575)   0.0 E45-E30 S56-S51
     17 1497  430 (S23W53 L= 51)   343     79        1 (1458/ 419-1541/ 444)  26.6 W50-W56 S24-S21
     18 1393  298 (S34W57 L= 56)  1094    142        1 (1340/ 250-1447/ 347) -82.9 W54-W61 S37-S29
  

    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