Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1017 1863 (N51E25 L=293)   785     82        3 ( 979/1839-1065/1889)  -5.2 E30-E19 N48-N53
      2  745 1702 (N31E34 L=283)  1325    321        2 ( 693/1552- 807/1823) -67.8 E41-E30 N19-N41
      3  956 1689 (N34E19 L=299)   370    104        1 ( 938/1639- 975/1747) -53.1 E22-E16 N30-N38
      4 1493 1514 (N33W22 L=340)   347     78        0 (1454/1497-1531/1532)  33.7 W19-W25 N30-N34
      5 1416 1391 (N23W18 L=335)   268     64        1 (1389/1356-1446/1413) -36.9 W15-W19 N20-N23
      6  608  938 (S18E24 L=293)   680    184        2 ( 553/ 868- 667/ 988)  52.1 E27-E20 S23-S14
      7  423  939 (S21E38 L=280)   239     47        0 ( 393/ 920- 450/ 958)  31.0 E40-E35 S22-S19
      8  413  581 (S44E42 L=276)  1240    161        4 ( 363/ 536- 466/ 637) -35.5 E46-E35 S47-S40
      9 1719  483 (S22W68 L= 25)  2483    312        4 (1602/ 412-1817/ 549)  40.9 W58-W78 S29-S15
     10 1095  217 (S57W36 L=353)   880    142        1 (1028/ 182-1148/ 246)  -3.0 W26-W45 S59-S54
  

    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