Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1500 1616 (N25W49 L= 83)  3211    553        7 (1271/1493-1732/1688) -32.3 W32-W63 N11-N31
      2 1713 1615 (N18W71 L=105)  3084    223        4 (1615/1543-1799/1680) -47.5 W63-W81 N12-N24
      3  674 1572 (N45E10 L= 24)   234     37        0 ( 649/1558- 703/1588) -26.6 E12-E07 N43-N46
      4 1390 1506 (N21W35 L= 69)   245     77        0 (1352/1494-1433/1516) -31.0 W32-W38 N19-N22
      5 1731 1214 (S03W49 L= 83)   300     35        1 (1711/1186-1752/1234)  33.7 W47-W51 S04-S01
      6 1731 1118 (S09W46 L= 80)  1443    355        3 (1634/1049-1794/1183) -39.3 W39-W51 S12-S03
      7  492 1102 (N19E31 L=  4)   816    276        1 ( 406/1025- 614/1176) -45.0 E36-E23 N12-N25
      8  916  874 (S03E09 L= 25)  1665    500        5 ( 763/ 687-1060/1002) -59.9 E16-E05 S16-N05
      9 1349  908 (S11W16 L= 50)   764    225        3 (1290/ 800-1392/ 992)  48.4 W10-W18 S15-S06
     10 1226  900 (S08W08 L= 43)   381    129        1 (1200/ 835-1246/ 962)  60.3 W05-W10 S11-S04
     11 1275  929 (S08W12 L= 46)   220     44        1 (1256/ 905-1297/ 956)  45.0 W10-W13 S09-S06
     12  390  906 (N10E42 L=352)   624    193        2 ( 314/ 874- 479/ 953) -26.6 E47-E35 N07-N14
     13 1469  678 (S27W19 L= 53)  1151    252        5 (1368/ 638-1573/ 728) -18.4 W12-W27 S29-S22
     14  710  659 (S10E27 L=  7)   299     72        2 ( 681/ 627- 737/ 703)  78.7 E27-E25 S12-S07
     15 1213  502 (S32E02 L= 33)  2025    458        3 (1135/ 318-1271/ 696)  62.5 E13-W03 S43-S18
     16  564  600 (S10E38 L=356)   571    124        4 ( 516/ 561- 614/ 643) -63.4 E40-E35 S13-S07
     17  653  587 (S13E33 L=  1)   379    111        1 ( 604/ 558- 695/ 621) -51.3 E34-E30 S15-S10
     18  535  418 (S21E49 L=345)  2771    477        7 ( 391/ 315- 656/ 563)  13.1 E77-E34 S25-S14
     19 1067  271 (S45E21 L= 13)  1717    242       10 (1013/ 194-1126/ 384)  39.1 E30-E14 S49-S36
  

    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