Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1114 1828 (N59W35 L= 57)  1503    254        3 ( 996/1783-1224/1866)   0.0 W18-W52 N56-N60
      2 1430 1596 (N34W43 L= 64)   519    185        0 (1341/1555-1505/1620) -38.7 W36-W47 N29-N37
      3 1362 1492 (N29W32 L= 54)   464    151        1 (1302/1472-1437/1517) -29.1 W28-W37 N26-N31
      4  463 1344 (N33E37 L=345)  1913    404        1 ( 405/1169- 527/1512)  72.3 E40-E32 N22-N44
      5 1829 1292 (N07W65 L= 86)   974    130        1 (1805/1240-1854/1365)  45.0 W60-W69 N04-N11
      6  259 1240 (N28E56 L=325)  2846    297        8 ( 167/1133- 351/1352) -25.8 E77-E45 N21-N35
      7  722 1212 (N22E16 L=  5)   294     89        0 ( 682/1197- 758/1237)   0.0 E19-E13 N20-N23
      8 1110 1146 (N12W07 L= 28)   243     96        0 (1068/1129-1159/1167) -31.0 W04-W09 N10-N14
      9 1569 1051 (S01W34 L= 56)   538     87        0 (1512/1033-1620/1070)  -7.1 W30-W38 S02-N00
     10 1701  864 (S15W43 L= 64)  1329    198        2 (1616/ 814-1775/ 928)  11.3 W35-W50 S17-S12
     11  868  622 (S16E16 L=  6)  1009    340        0 ( 747/ 555- 951/ 735)  32.7 E24-E09 S19-S09
     12  551  639 (S11E37 L=345)  1128    189        2 ( 531/ 538- 588/ 718)  61.4 E40-E32 S17-S05
     13  700  643 (S12E26 L=355)   326     60        2 ( 662/ 629- 741/ 660) -11.3 E28-E23 S13-S10
     14  277  535 (S16E69 L=313)   739    122        2 ( 249/ 503- 292/ 589)  33.7 E75-E63 S18-S12
     15  321  547 (S15E61 L=320)   517     51        1 ( 306/ 513- 336/ 581)  59.0 E64-E58 S17-S12
     16 1548  398 (S46W33 L= 55)  1610    281        4 (1457/ 339-1682/ 495)  15.3 W24-W46 S49-S40
     17 1347  315 (S48W13 L= 35)   424     62        1 (1324/ 290-1375/ 340) -59.0 W11-W16 S50-S45
     18 1280  251 (S53W06 L= 28)  1068    114        3 (1230/ 220-1328/ 292) -11.3 W01-W11 S57-S48
     19 1276  189 (S61W06 L= 27)   647     70        1 (1235/ 173-1311/ 206) -26.6 W00-W10 S63-S57
  

    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