Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1204 1842 (N54E22 L=292)   564     73        0 (1174/1821-1234/1862)  -8.1 E25-E17 N50-N55
      2  628 1800 (N30E60 L=255)  1352    153        1 ( 588/1738- 684/1855) -54.0 E63-E55 N24-N35
      3  854 1767 (N34E38 L=276)  1162    226        2 ( 829/1664- 887/1847) -35.2 E46-E29 N27-N39
      4  864 1446 (N13E21 L=293)  1148    294        4 ( 829/1316- 897/1589) -62.1 E27-E17 N04-N21
      5  387 1399 (S01E51 L=263)   302     43        0 ( 365/1375- 409/1418)  63.4 E52-E50 S02-N00
      6  218 1097 (S22E59 L=255)   846    226        1 ( 185/1003- 250/1201) -67.4 E64-E54 S27-S15
      7 1076  625 (S27W16 L=330)   218     80        0 (1047/ 589-1093/ 657)  90.0 W14-W16 S29-S24
      8 1413  619 (S16W36 L=350)   433     56        0 (1389/ 597-1449/ 652)   0.0 W33-W38 S18-S14
      9 1233  522 (S27W30 L=344)   875    173        1 (1155/ 504-1317/ 547)  41.2 W24-W33 S30-S23
     10 1491  458 (S22W50 L=  5)  1031    123        3 (1425/ 423-1546/ 504)  50.2 W46-W53 S25-S18
     11  680  308 (S64W02 L=317)  1392    145        4 ( 623/ 275- 745/ 363)  14.9 E04-W11 S67-S59
     12  824  320 (S56W15 L=330)   495     63        1 ( 805/ 286- 846/ 355) -50.2 W12-W18 S58-S52
  

    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