Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  864 1807 (N49E37 L=254)   542     77        0 ( 834/1778- 906/1825)  68.2 E39-E33 N46-N51
      2  656 1717 (N37E47 L=245)   899    138        1 ( 627/1659- 696/1785) -27.9 E57-E39 N33-N42
      3  972 1684 (N41E19 L=272)   364     54        0 ( 946/1649- 993/1710) -38.7 E21-E16 N38-N43
      4  450 1638 (N28E60 L=231)   463     55        2 ( 427/1618- 473/1668) -18.4 E65-E56 N26-N29
      5 1557 1559 (N44W34 L=325)   559     84        1 (1515/1535-1597/1579)   6.3 W29-W38 N42-N46
      6 1381 1457 (N33W17 L=308)   437    144        1 (1326/1410-1437/1503) -29.1 W12-W21 N30-N35
      7 1084 1363 (N21E02 L=289)   431    169        0 (1044/1308-1143/1419) -35.5 E05-W02 N18-N23
      8  282 1211 (S02E55 L=236)   813    159        1 ( 233/1168- 333/1269) -18.4 E61-E49 S03-N01
      9 1782  647 (S10W65 L=356)   610     80        0 (1736/ 626-1820/ 672)  21.8 W59-W69 S11-S07
     10  937  562 (S31W03 L=294)   329     77        1 ( 899/ 544- 976/ 580)  31.0 W00-W05 S32-S28
     11  984  373 (S44W12 L=303)   423     99        0 ( 946/ 335-1019/ 406)  74.1 W09-W14 S47-S40
     12 1546  338 (S34W63 L=354)  1113    124        2 (1500/ 317-1590/ 373)  10.3 W57-W70 S36-S31
  

    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