Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  838 1657 (N32E21 L= 94)   356     47        1 ( 817/1636- 871/1686)   0.0 E22-E18 N30-N34
      2 1252 1563 (N30W09 L=124)   513    170        0 (1208/1492-1293/1637) -52.1 W05-W12 N25-N35
      3 1052 1437 (N19E03 L=112)   404    135        0 (1037/1370-1080/1500) -63.4 E04-W00 N15-N23
      4 1485 1464 (N27W26 L=141)   286     79        0 (1453/1442-1525/1489) -26.6 W23-W29 N25-N28
      5  877 1325 (N11E12 L=103)   208     85        0 ( 848/1293- 906/1358)  59.0 E13-E10 N08-N13
      6 1762 1221 (N15W48 L=163)   683     88        2 (1740/1177-1804/1261)  63.4 W45-W52 N12-N18
      7 1698 1111 (N08W42 L=157)  1775    414        2 (1561/ 948-1777/1248) -34.4 W31-W50 N00-N14
      8 1843 1206 (N16W56 L=171)   414     40        0 (1829/1177-1859/1237) -90.0 W54-W58 N14-N18
      9 1841 1108 (N10W56 L=171)   339     61        1 (1813/1086-1866/1133)  18.4 W52-W58 N08-N11
     10  306 1021 (S11E48 L= 67)   554     78        0 ( 285/ 972- 324/1067) -78.7 E49-E46 S13-S08
     11  191 1023 (S11E60 L= 55)   720     71        1 ( 156/ 989- 221/1055) -23.2 E64-E56 S13-S09
     12 1243 1023 (S03W13 L=128)   505    177        0 (1155/ 963-1292/1054) -12.5 W07-W16 S06-S01
     13  232  854 (S21E56 L= 59)  3315    447        3 ( 121/ 762- 341/ 953)  27.6 E70-E44 S26-S15
     14  517  871 (S19E31 L= 84)   937    315        2 ( 445/ 813- 583/ 952) -49.4 E36-E26 S22-S14
     15 1801  717 (S13W60 L=175)  6772      0       27 (1568/ 552-1937/ 883)  25.2 W41-W76 S23-S02
     16 1194  844 (S14W12 L=127)   195     53        1 (1165/ 826-1226/ 865) -14.0 W10-W14 S15-S13
     17 1482  655 (S22W35 L=150)   622    152        0 (1458/ 592-1526/ 733) -41.2 W31-W39 S25-S17
     18  440  549 (S40E42 L= 73)  1596    209        4 ( 342/ 516- 550/ 587)   7.1 E54-E30 S42-S37
  

    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