Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  989 1680 (N49W00 L=169)   361     52        1 ( 968/1648-1016/1707)  68.2 E01-W02 N46-N51
      2 1140 1664 (N47W14 L=183)   664    164        0 (1058/1639-1218/1684)  -8.1 W06-W20 N44-N48
      3  384 1568 (N41E59 L=110)   451     87        0 ( 348/1556- 433/1581)   5.2 E65-E52 N39-N41
      4 1359 1500 (N33W27 L=196)   266     61        0 (1329/1477-1388/1522) -26.6 W25-W29 N31-N34
      5  388 1343 (N25E47 L=122)   461    157        0 ( 340/1289- 424/1402)  58.0 E50-E43 N21-N29
      6  696 1290 (N22E21 L=148)   586    159        1 ( 628/1245- 745/1325) -20.6 E25-E17 N18-N23
      7 1406 1266 (N17W26 L=195)   399    144        0 (1355/1217-1447/1322) -54.5 W23-W28 N14-N21
      8  517 1216 (N17E34 L=135)   254     63        1 ( 492/1190- 538/1248)  45.0 E35-E32 N15-N19
      9 1201  953 (S01W11 L=179)   551    266        1 (1095/ 881-1303/ 998) -26.6 W04-W16 S05-N02
     10  376  708 (S14E47 L=122)  1442    378        4 ( 338/ 554- 420/ 864)  63.4 E53-E40 S24-S04
     11 1700  778 (S15W47 L=216)   424    109        2 (1650/ 742-1726/ 819)  38.7 W42-W49 S16-S12
     12  974  664 (S18E05 L=164)  1405    328        6 ( 834/ 588-1060/ 728) -37.6 E13-W00 S23-S13
     13 1665  607 (S26W47 L=216)   821    202        0 (1606/ 543-1718/ 685)  52.1 W43-W50 S30-S21
     14 1292  599 (S24W16 L=185)  1270    270        3 (1203/ 502-1397/ 676)  37.6 W10-W23 S30-S19
     15  512  312 (S44E56 L=113)   585     53        2 ( 488/ 295- 539/ 338) -71.6 E60-E52 S46-S41
     16 1115  234 (S54W05 L=173)   454     74        0 (1098/ 209-1139/ 268) -63.4 W02-W07 S56-S50
     17 1225  211 (S58W18 L=187)   838     87        1 (1192/ 179-1255/ 239) -45.0 W13-W22 S62-S53
     18 1302  189 (S62W31 L=200)  1117     73        2 (1283/ 147-1321/ 232) -31.0 W24-W49 S70-S55
  

    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