Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1552 1667 (N39W45 L= 77)   391     42        0 (1528/1646-1573/1690)  23.2 W40-W47 N37-N40
      2 1331 1598 (N33W22 L= 54)   273     84        0 (1287/1581-1370/1613) -18.4 W18-W24 N31-N34
      3 1820 1460 (N25W65 L= 97)   900     75        1 (1780/1420-1858/1502) -35.5 W61-W70 N22-N27
      4  502 1268 (N12E33 L=359)   723    135        2 ( 478/1211- 547/1338) -63.4 E35-E29 N08-N16
      5  328 1233 (N10E47 L=345)   637    143        1 ( 298/1166- 372/1289)  54.5 E48-E43 N06-N13
      6  390 1220 (N09E41 L=351)   293     32        1 ( 371/1199- 409/1242)  71.6 E42-E39 N07-N10
      7  720 1082 (N00E18 L= 14)   825    252        2 ( 609/1038- 828/1118)   4.1 E25-E11 S02-N02
      8 1043  888 (S12W01 L= 33)  1560    501        3 ( 989/ 665-1085/1099) -87.8 E02-W03 S25-N01
      9 1789  922 (S08W53 L= 85)   885    180        3 (1758/ 859-1818/1008)  82.9 W49-W55 S11-S03
     10  637  893 (S11E24 L=  8)   740    140        5 ( 575/ 855- 698/ 943)  12.5 E28-E19 S13-S07
     11 1724  903 (S10W47 L= 79)   428     35        1 (1702/ 869-1747/ 932) -76.0 W44-W49 S11-S07
     12 1767  794 (S16W53 L= 85)  1559    321        2 (1676/ 723-1872/ 867)  26.6 W45-W63 S20-S11
     13 1576  763 (S18W37 L= 69)   795    243        0 (1518/ 666-1606/ 866)  76.0 W33-W39 S24-S12
     14  603  740 (S20E28 L=  5)   371    142        1 ( 565/ 702- 667/ 787) -26.6 E29-E23 S22-S17
     15  415  654 (S25E44 L=348)   979    179        4 ( 334/ 595- 480/ 714)  21.8 E52-E37 S28-S21
     16 1417  668 (S25W26 L= 59)   233     39        1 (1396/ 651-1446/ 689) -14.0 W24-W28 S25-S23
     17  927  592 (S30E07 L= 25)   557    145        0 ( 858/ 554- 989/ 634)  24.0 E11-E02 S32-S27
     18  503  417 (S41E46 L=347)  3550    361        7 ( 328/ 367- 673/ 477)   5.4 E80-E27 S43-S37
     19 1554  364 (S45W51 L= 83)  3317    405        8 (1381/ 260-1707/ 447)  30.0 W39-W65 S54-S38
  

    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