Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  503 1673 (N20E56 L=314)  2166    273        6 ( 431/1625- 572/1735)  63.4 E61-E52 N15-N25
      2 1337 1591 (N38W03 L= 13)   349     54        0 (1302/1566-1362/1611) -14.0 W00-W04 N36-N39
      3  766 1512 (N16E29 L=341)   456    161        0 ( 743/1459- 808/1584) -36.9 E33-E24 N13-N20
      4 1430 1464 (N33W13 L= 23)  1581    240        4 (1370/1351-1475/1558) -50.2 W07-W18 N26-N38
      5  492 1452 (N06E45 L=325)   292     47        0 ( 466/1424- 519/1475)  63.4 E46-E43 N03-N07
      6  573 1406 (N05E37 L=332)   381    137        1 ( 507/1360- 614/1465)  69.4 E40-E35 N01-N09
      7 1800 1159 (N26W48 L= 57)   704    148        0 (1778/1090-1816/1231) -76.0 W44-W49 N21-N30
      8 1747 1024 (N16W45 L= 55)   401     82        3 (1724/1000-1783/1050)  11.3 W43-W48 N14-N17
      9 1905  921 (N16W65 L= 75)  2316    258        4 (1884/ 827-1926/1020) -63.4 W60-W73 N10-N22
     10  924  930 (S13E03 L=  7)   195    103        0 ( 865/ 898- 973/ 963)  56.3 E05-E01 S15-S09
     11  258  889 (S32E50 L=320)   744    213        1 ( 229/ 809- 319/ 937)  59.0 E52-E44 S37-S27
     12 1605  870 (N03W39 L= 49)   463    142        2 (1537/ 837-1657/ 915)  35.5 W34-W42 N00-N06
     13  308  763 (S39E44 L=326)   617    135        2 ( 274/ 732- 354/ 802) -26.6 E47-E39 S41-S36
     14 1433  692 (S12W34 L= 43)   550    195        1 (1345/ 621-1519/ 732)  55.0 W30-W38 S17-S07
     15 1151  706 (S19W16 L= 26)   334    154        0 (1074/ 684-1227/ 727)  16.7 W11-W21 S21-S17
     16 1197  571 (S25W24 L= 34)   576    161        1 (1115/ 558-1276/ 593)  29.1 W19-W28 S27-S22
  

    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