Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1082 1762 (N45E06 L=261)   859    107        0 (1023/1737-1146/1787)  -5.2 E11-W00 N43-N46
      2 1534 1629 (N41W34 L=301)   426     77        0 (1489/1606-1570/1646)   0.0 W29-W37 N39-N42
      3  786 1469 (N19E20 L=247)   978    254        0 ( 706/1369- 877/1549) -40.2 E27-E12 N13-N24
      4 1247 1466 (N24W10 L=277)   246     75        0 (1217/1441-1282/1494) -26.6 W07-W11 N23-N25
      5 1482 1430 (N25W26 L=293)   367     43        0 (1458/1399-1508/1461) -45.0 W24-W28 N23-N27
      6 1518 1264 (N15W29 L=296)   216     38        1 (1503/1232-1532/1289)  90.0 W27-W29 N12-N16
      7 1680 1224 (N15W41 L=308)   574    158        0 (1641/1155-1717/1288) -49.4 W38-W44 N11-N18
      8  224  746 (S28E62 L=205)   974    133        0 ( 202/ 685- 258/ 808)  30.3 E68-E55 S32-S24
      9  980  644 (S29W01 L=269)  1136    363        0 ( 856/ 527-1102/ 743)  48.8 E06-W08 S37-S21
     10 1732  691 (S16W55 L=322)   573     58        0 (1710/ 658-1757/ 726)  63.4 W53-W57 S17-S13
     11 1791  678 (S15W62 L=329)   692    119        3 (1756/ 637-1828/ 704)   0.0 W58-W67 S17-S13
  

    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