Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  772 1734 (N58W09 L=345)  2946    510        7 ( 598/1590- 970/1873)   7.4 E14-W42 N52-N62
      2 1573 1681 (N25W65 L= 40)  3657    354       10 (1446/1605-1716/1756) -50.5 W56-W76 N15-N33
      3 1601 1577 (N19W58 L= 33)   681    155        0 (1528/1559-1678/1601) -32.5 W52-W64 N15-N22
      4 1338 1515 (N24W35 L= 10)   246     86        1 (1300/1493-1379/1542)  -9.5 W32-W38 N23-N25
      5  731 1315 (N30E10 L=326)   319     84        2 ( 685/1288- 768/1347)   8.1 E13-E06 N29-N31
      6  965 1221 (N18W02 L=338)   219     66        1 ( 942/1191- 990/1257)  45.0 W00-W03 N16-N20
      7 1820  837 (S29W50 L= 26)   600     74        1 (1788/ 820-1859/ 858) -29.1 W46-W55 S31-S26
      8 1413  752 (S20W15 L=351)   426    113        0 (1394/ 693-1427/ 813)  60.3 W13-W17 S23-S16
      9 1241  634 (S22W02 L=337)   283    125        0 (1210/ 585-1264/ 696)  50.2 E01-W04 S24-S18
     10 1765  611 (S43W45 L= 20)   474     45        1 (1735/ 590-1790/ 632)  -7.1 W40-W48 S44-S41
     11 1349  536 (S32W06 L=341)   279     54        0 (1331/ 507-1368/ 560)  45.0 W04-W07 S33-S30
     12 1337  399 (S42W01 L=337)   640    139        2 (1310/ 361-1384/ 470) -90.0 E01-W03 S45-S35
     13  912  394 (S28E26 L=309)   446    116        1 ( 872/ 362- 970/ 425) -66.8 E28-E24 S31-S24
     14  534  317 (S24E60 L=276)   629    110        1 ( 486/ 279- 573/ 346) -63.4 E62-E56 S27-S21
     15  686  283 (S30E49 L=286)   492     87        1 ( 663/ 252- 706/ 323)  45.0 E53-E45 S32-S26
     16  691  226 (S34E55 L=280)   534     82        3 ( 642/ 209- 730/ 237) -12.5 E59-E50 S36-S32
  

    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