Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  978 1690 (N36E23 L=135)  2501    316        7 ( 869/1565-1064/1784) -24.6 E35-E11 N29-N42
      2  543 1714 (N27E57 L=101)   505     72        0 ( 521/1674- 570/1748) -20.0 E62-E51 N25-N29
      3 1607 1599 (N50W32 L=190)  1531    185        1 (1542/1557-1697/1643)  -2.9 W24-W44 N47-N52
      4  422 1321 (N01E45 L=113)  1478    218        3 ( 385/1239- 475/1405) -45.0 E49-E39 S03-N05
      5 1291 1326 (N21W09 L=166)   528    173        0 (1224/1261-1341/1383) -29.1 W03-W12 N17-N22
      6  541  962 (S17E28 L=130)   261     47        0 ( 523/ 937- 565/ 984)  90.0 E29-E26 S18-S15
      7 1039  830 (S14W05 L=163)  1039    318        0 ( 903/ 759-1177/ 877)  45.0 W00-W12 S20-S08
      8 1494  790 (S05W34 L=191)   370    127        0 (1456/ 732-1515/ 838) -39.8 W29-W35 S08-S03
      9 1168  602 (S24W20 L=177)  1211    264        2 (1054/ 580-1303/ 642)  20.6 W12-W28 S26-S20
     10  454  567 (S45E32 L=125)   562     47        2 ( 419/ 542- 491/ 593)  20.6 E36-E28 S47-S42
     11 1535  367 (S28W59 L=217)   435     70        0 (1518/ 343-1566/ 397)  -7.1 W55-W63 S29-S25
  

    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