Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  783 1817 (N57E39 L=308)  1315    137        2 ( 740/1755- 832/1852) -54.0 E46-E33 N50-N61
      2  662 1822 (N56E58 L=290)   617     47        1 ( 634/1806- 690/1838)  -4.4 E65-E50 N54-N57
      3  968 1595 (N39E09 L=338)   254     48        0 ( 949/1570- 993/1610)  33.7 E10-E07 N37-N40
      4  531 1489 (N28E41 L=306)   699    188        1 ( 456/1454- 615/1537) -13.2 E50-E33 N26-N30
      5 1509 1409 (N29W33 L= 20)   673    207        1 (1451/1336-1577/1480) -48.4 W29-W37 N24-N34
      6  326 1392 (N20E57 L=291)  1385    167        3 ( 269/1321- 366/1455) -22.4 E67-E50 N15-N23
      7  805 1303 (N18E16 L=331)   635    272        2 ( 703/1213- 906/1363) -23.6 E24-E08 N12-N21
      8 1353 1074 (N07W20 L=  8)   423     87        2 (1324/1036-1379/1111)  59.0 W18-W22 N04-N09
      9 1270  699 (S17W18 L=  5)   220     72        0 (1249/ 659-1290/ 731)  78.7 W17-W19 S19-S14
     10 1265  297 (S47W30 L= 17)   473     37        0 (1242/ 271-1286/ 318) -26.6 W26-W32 S49-S44
  

    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