Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1496 1563 (N21W48 L=208)  1067    192        3 (1419/1509-1587/1595) -10.6 W39-W55 N17-N22
      2 1079 1527 (N32W20 L=180)   411     74        0 (1039/1503-1120/1548)  -9.5 W16-W22 N30-N33
      3 1715 1418 (N06W57 L=218)   380     84        2 (1686/1396-1750/1445) -53.1 W55-W60 N04-N08
      4  621 1376 (N37E16 L=144)   235     50        0 ( 589/1365- 655/1387) -11.3 E18-E13 N35-N38
      5 1718 1312 (N01W52 L=213)   844    122        0 (1690/1245-1754/1369) -90.0 W50-W54 S03-N04
      6 1212  708 (S17W02 L=162)   473    201        0 (1130/ 661-1268/ 777)  18.4 E03-W06 S18-S14
      7  795  651 (S09E23 L=137)   856    297        3 ( 714/ 530- 866/ 769)  86.6 E25-E21 S17-N00
      8 1742  650 (S39W41 L=201)  1297    146        0 (1708/ 578-1795/ 695)  36.9 W37-W47 S43-S35
      9 1580  571 (S38W24 L=184)   277     43        1 (1563/ 545-1597/ 596)  45.0 W22-W25 S39-S36
     10 1284  535 (S30W01 L=162)   420    119        0 (1230/ 518-1336/ 563)  -9.5 E01-W05 S31-S27
     11  740  482 (S18E33 L=127)   576    184        0 ( 684/ 447- 777/ 547)  29.7 E36-E28 S19-S14
     12  946  327 (S34E28 L=132)   790    127        2 ( 917/ 264- 973/ 378)  37.9 E33-E24 S38-S29
  

    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