Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  822 1864 (N67W21 L= 26)   997    109        2 ( 795/1820- 853/1910)  18.4 W12-W36 N62-N71
      2  910 1880 (N64W33 L= 37)   446     49        0 ( 897/1851- 923/1909)  30.3 W27-W40 N60-N67
      3 1165 1753 (N43W36 L= 41)   766    132        1 (1130/1708-1196/1806)  28.6 W30-W41 N40-N47
      4  929 1748 (N51W18 L= 23)   488     59        0 ( 888/1731- 967/1763) -18.4 W14-W21 N49-N53
      5 1454 1642 (N28W49 L= 53)   678    118        1 (1403/1615-1514/1680) -60.3 W45-W52 N24-N31
      6 1288 1567 (N28W32 L= 37)   852    329        1 (1168/1523-1429/1642) -47.3 W27-W40 N22-N35
      7 1574 1429 (N12W45 L= 50)   677    148        2 (1535/1350-1608/1493)  66.0 W43-W47 N07-N16
      8  920 1427 (N29W04 L=  9)   453    163        0 ( 845/1381- 976/1447) -49.4 W00-W06 N25-N32
      9 1195 1425 (N22W21 L= 25)   415    227        0 (1082/1413-1295/1438) -24.8 W13-W26 N18-N25
     10 1133 1362 (N20W15 L= 19)   322    129        0 (1083/1339-1194/1388) -45.0 W12-W17 N17-N22
     11  349 1124 (N24E42 L=322)   484     80        0 ( 318/1085- 388/1160) -50.2 E45-E39 N21-N27
     12  313 1001 (N18E47 L=317)   635     67        1 ( 282/ 971- 347/1034) -45.0 E50-E44 N15-N20
     13  209  876 (N12E60 L=304)   390     77        0 ( 193/ 834- 224/ 912)  45.0 E63-E58 N10-N14
     14  833  664 (S12E20 L=345)   383    106        1 ( 809/ 604- 865/ 705)  39.8 E22-E16 S15-S09
     15 1709  649 (S37W40 L= 44)   448     23        0 (1688/ 625-1728/ 673) -63.4 W37-W42 S38-S34
     16 1294  446 (S38W02 L=  7)  1678    357        5 (1173/ 325-1403/ 557)  22.2 E10-W12 S44-S32
     17 1458  272 (S61W14 L= 19)   711     59        2 (1418/ 251-1499/ 295) -14.0 W08-W20 S63-S58
     18 1322  186 (S65E07 L=358)   517     58        0 (1293/ 168-1350/ 205) -14.0 E10-E01 S67-S61
  

    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