Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  407 1643 (N23E67 L=231)  1567    118        3 ( 386/1583- 442/1697) -24.0 E80-E60 N18-N26
      2  618 1645 (N27E45 L=252)   351     89        0 ( 576/1623- 654/1664)   0.0 E49-E41 N25-N28
      3 1227 1436 (N27W04 L=301)   915    359        1 (1085/1332-1351/1495) -11.9 E06-W13 N23-N29
      4 1398 1432 (N31W16 L=313)   398    155        0 (1347/1371-1448/1491) -29.1 W11-W20 N28-N34
      5 1893 1170 (N27W67 L=  5)  1785    207        0 (1876/1077-1907/1272)  77.0 W64-W70 N20-N33
      6  480 1057 (S11E35 L=263)   390     77        0 ( 440/1037- 514/1078)  31.0 E37-E32 S12-S09
      7 1791  985 (N12W54 L=351)   666     93        0 (1767/ 933-1811/1036)  90.0 W51-W55 N09-N15
      8 1472  548 (S20W41 L=339)  1923    362        2 (1394/ 428-1588/ 636)  79.4 W38-W46 S28-S12
      9  491  541 (S44E31 L=267)   373     47        0 ( 476/ 510- 507/ 575) -90.0 E32-E28 S46-S41
     10  933  516 (S35W06 L=304)   622    174        1 ( 852/ 488-1022/ 547)  28.6 W00-W11 S37-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