Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1311 1658 (N42W34 L=287)   276     81        0 (1275/1646-1355/1670) -26.6 W31-W37 N40-N43
      2 1523 1539 (N30W45 L=298)   437     90        2 (1484/1503-1551/1559) -53.1 W42-W47 N27-N31
      3  996 1335 (N25W02 L=255)   254     65        1 ( 969/1307-1030/1362) -53.1 W00-W04 N23-N27
      4 1814 1308 (N11W64 L=317)  1344    225        4 (1766/1259-1859/1353)  14.0 W57-W70 N08-N13
      5  539 1149 (N18E31 L=223)  1794    496        6 ( 375/1048- 706/1260) -33.7 E43-E19 N10-N26
      6 1846 1225 (N05W65 L=318)   526     81        3 (1829/1177-1864/1257)  78.7 W63-W67 N02-N07
      7 1310 1168 (N11W20 L=273)   396    140        2 (1253/1112-1351/1222)  32.0 W15-W23 N08-N13
      8 1172 1130 (N10W10 L=263)   244    115        0 (1110/1110-1228/1155) -29.7 W06-W13 N08-N12
      9 1844 1036 (S06W60 L=313)  1281    200        0 (1820/ 950-1865/1100)  71.6 W57-W63 S11-S02
     10  374 1002 (N10E44 L=209)   635    181        0 ( 324/ 925- 425/1055) -48.8 E48-E40 N05-N13
     11 1887  979 (S10W65 L=318)   556    111        2 (1868/ 945-1909/1010) -51.3 W62-W68 S13-S08
     12 1600  738 (S20W35 L=288)  2366    449        9 (1455/ 596-1696/ 896)  41.6 W24-W43 S27-S11
     13  901  774 (S08E11 L=243)   197     45        1 ( 881/ 752- 929/ 793)  18.4 E11-E08 S09-S06
     14  972  695 (S14E07 L=246)   844    328        1 ( 843/ 621-1069/ 786)  23.6 E16-W00 S17-S08
     15  791  708 (S11E18 L=235)   344     83        1 ( 755/ 677- 815/ 752)  63.4 E20-E16 S12-S08
     16  661  696 (S11E27 L=226)   506    130        5 ( 616/ 653- 717/ 749) -56.3 E29-E23 S13-S06
     17 1269  593 (S24W11 L=264)   226     52        1 (1253/ 570-1289/ 621) -90.0 W09-W12 S26-S22
     18  538  583 (S17E39 L=215)   656    144        2 ( 488/ 546- 584/ 619) -51.3 E41-E35 S19-S14
  

    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