Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  748 1545 (N35E24 L=134)   352     77        0 ( 709/1510- 776/1568) -23.2 E28-E21 N32-N36
      2 1585 1517 (N37W44 L=203)  1204    175        2 (1507/1476-1658/1544)  -3.6 W36-W52 N34-N38
      3 1326 1450 (N31W20 L=178)   294     92        1 (1296/1426-1361/1476) -36.9 W17-W22 N29-N32
      4  362 1448 (N25E55 L=103)   437     64        0 ( 332/1426- 393/1466)  31.0 E58-E52 N23-N26
      5 1800 1322 (N24W62 L=221)  1883    262        5 (1703/1250-1887/1392) -22.2 W52-W74 N19-N28
      6 1407 1365 (N26W25 L=184)   403     82        0 (1371/1328-1435/1391) -11.3 W22-W27 N23-N27
      7 1750  858 (S05W53 L=211)   447      7        0 (1727/ 839-1772/ 879) -11.3 W50-W55 S06-S03
      8  315  660 (S25E54 L=105)  2343    348        4 ( 240/ 533- 399/ 766) -56.3 E64-E48 S33-S18
      9 1647  687 (S16W46 L=205)   318     67        1 (1630/ 650-1665/ 720) -76.0 W45-W47 S18-S14
     10 1764  603 (S22W63 L=222)  1301    158        2 (1738/ 541-1790/ 681) -63.4 W59-W71 S26-S16
     11  597  511 (S34E30 L=129)   474     64        0 ( 560/ 492- 638/ 535)   0.0 E33-E26 S34-S31
     12  837  515 (S32E11 L=148)   231     43        0 ( 807/ 496- 865/ 531) -26.6 E12-E08 S32-S30
  

    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