Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1209 1839 (N53E21 L=292)   593     74        0 (1177/1814-1240/1861) -12.5 E25-E16 N50-N55
      2  635 1799 (N30E59 L=255)  1107    136        3 ( 594/1744- 690/1855) -55.0 E63-E54 N25-N35
      3  859 1764 (N33E38 L=276)  1177    226        2 ( 834/1663- 892/1848) -37.4 E46-E29 N27-N40
      4  870 1442 (N13E21 L=293)  1251    295        2 ( 830/1311- 899/1587) -64.8 E26-E16 N04-N21
      5  392 1396 (S01E51 L=263)   320     53        0 ( 366/1373- 415/1418)  33.7 E52-E49 S02-N00
      6  234 1025 (S26E56 L=258)   342     43        1 ( 218/ 998- 252/1051)  56.3 E57-E54 S27-S24
      7 1415  616 (S16W36 L=350)   482     65        1 (1379/ 593-1454/ 649)  36.9 W34-W39 S18-S14
      8 1236  519 (S27W30 L=344)   898    175        2 (1149/ 499-1319/ 542)  37.9 W24-W34 S30-S23
      9 1491  452 (S22W51 L=  4)  1141    154        1 (1411/ 417-1550/ 502)  45.0 W46-W54 S26-S18
     10  515  373 (S65E23 L=290)   394     34        1 ( 497/ 350- 534/ 391)  45.0 E26-E19 S67-S62
     11  820  315 (S56W15 L=329)   666     90        1 ( 785/ 283- 850/ 363)  81.9 W11-W18 S60-S52
     12  677  307 (S64W02 L=316)  1454    115        2 ( 625/ 276- 740/ 359)  33.7 E04-W09 S68-S59
  

    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