Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1445 1515 (N34W33 L= 89)   486     89        0 (1421/1474-1467/1560) -90.0 W31-W34 N31-N37
      2 1661 1496 (N32W51 L=107)   843    109        2 (1621/1460-1698/1538) -51.3 W48-W54 N29-N34
      3  316 1484 (N32E57 L=359)  1064    236        3 ( 284/1443- 364/1533) -28.6 E62-E51 N29-N35
      4  574 1392 (N27E32 L= 24)   261     69        0 ( 542/1380- 609/1405) -11.3 E34-E29 N26-N27
      5  172 1223 (N15E68 L=348)   635     87        0 ( 156/1194- 188/1255) -36.9 E70-E64 N13-N16
      6  752 1179 (N14E17 L= 39)   236     63        1 ( 718/1152- 784/1207) -45.0 E19-E15 N11-N15
      7 1149  969 (N00W08 L= 64)   290     64        1 (1108/ 945-1190/ 994)  11.3 W05-W10 S01-N01
      8 1886  955 (S04W67 L=123)   573    121        0 (1863/ 928-1914/ 977)  20.6 W64-W72 S05-S02
      9 1850  832 (S11W63 L=119)   790    161        0 (1836/ 796-1879/ 894) -63.4 W60-W67 S13-S07
     10  562  747 (S13E31 L= 25)   445    112        0 ( 503/ 719- 607/ 780) -23.2 E34-E27 S15-S11
     11  378  618 (S22E48 L=  8)  1575    388        1 ( 293/ 530- 448/ 703)  35.8 E60-E41 S29-S16
     12 1040  571 (S25W00 L= 56)   241     46        0 (1006/ 558-1068/ 585) -18.4 E01-W02 S26-S24
  

    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